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OPTIMAL REGULARITY AND EXPONENTIAL STABILITY FOR THE 
BLACKSTOCK-CRIGHTON EQUATION IN Lp-SPACES WITH 
DIRICHLET AND NEUMANN BOUNDARY CONDITIONS 

RAINER BRUNNHUBER AND STEFAN MEYER 


Abstract. The Blackstock-Crighton equation models nonlinear acoustic wave propagation 
in thermo-viscous fluids. In the present work we investigate the associated inhomogeneous 
Dirichlet and Neumann boundary value problems in a bounded domain and prove long-time 
well-posedness and exponential stability for sufficiently small data. The solution depends 
analytically on the data. In the Dirichlet case, the solution decays to zero and the same 
holds for Neumann conditions if the data have zero mean. 

We choose an optimal Lp-setting, where the regularity of the initial and boundary data 
are necessary and sufficient for existence, uniqueness and regularity of the solution. The 
linearized model with homogeneous boundary conditions is represented as an abstract evolu¬ 
tion equation for which we show maximal Lp-regularity. In order to eliminate inhomogeneous 
boundary conditions, we establish a general higher regularity result for the heat equation. 
We conclude that the linearized model induces a topological linear isomorphism and then 
solve the nonlinear problem by means of the implicit function theorem. 
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1. Introduction 

An acoustic wave propagates through a medium as a local pressure change. Nonlinear 
effects typically occur in case of acoustic waves of high amplitude which are used for several 
medical and industrial purposes such as lithotripsy, thermotherapy, ultrasound cleaning or 
welding and sonochemistry. Research on mathematical aspects of nonlinear acoustic wave 
propagation is therefore not only interesting from a mathematicians point of view. In fact, 
in case of medical applications, enhancement of the mathematical understanding of the un¬ 
derlying models should lead to a considerable reduction of complication risks. 
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The present work aims to provide a mathematical analysis of the Blackstock-Crighton- 
Kuznetsov equation 

(1.1) (oA - dt) {utt - c^Au - bAut) = + |Vu|^ 

and the Blackstock-Crighton-Westervelt equation 

(1.2) (oA - dt) {utt - c^Au - bAut) = (^ (l + {ut)‘^)tt 

for the acoustic velocity potential u, where c is the speed of sound, b is the diffusivity of 
sound and a is the heat conductivity of the fluid. Note that a = uPr, where u is the is 
kinematic viscosity and Pr denotes the Prandtl number. Alternatively, (11.11) and (|1.2I) can be 
expressed in terms of the acoustic pressure p via the pressure density relation put = p, where 
p denotes the mass density. The quantity B/A\s known as the parameter of nonlinearity 
and is proportional to the ratio of the coefficients of the quadratic and linear terms in the 
Taylor series expansion of the variations of the pressure in a medium in terms of variations 
of the density. Note that (jl.2|) is obtained from (jl.lj) by neglecting local nonlinear effects in 
the sense that the expression c^|Vup — {ut)"^ is sufficiently small. For a detailed introduction 
to the theory and applications of nonlinear acoustics we refer to [HB98| . 

Equations (II.ip and (II.2p result from two evolution equations of fourth order governing 
finite-amplitude sound in thermoviscous relaxing fluids, namely 

(1.3) -c^aA^u + {a + b) Auu + {c^Au - uut) = (|Vu|j -b ^utAu)^, 

(1.4) (aA - dt) {utt - c^Au) = (|Vu|^ -b ^utAu)^, 

which have been derived by Blackstock |Bla63| from the basic equations describing the gen¬ 
eral motion of thermally relaxing, viscous fluids (continuity equation, momentum equation, 
entropy equation and an arbitrary equation of state) and also appear as equations (11) and 
(13) in Crighton’s work |Cri79j on nonlinear acoustic models. We replace Au in the last 
term of and (jl.4p by -^Utt, which can be justified by the main part of the differential 
operator corresponding to the wave equation uu — cPAu = 0. Moreover, in (II.4p we consider 
potential diffusivity as in (11.31) . Therewith, we arrive at equation dni for which in [Brul5] 
the name Blackstock-Crighton-Kuznetsov equation has been introduced. For a more rigorous 
derivation of (11.11) we refer to Section 2 in |Bru IS]. 

While (|l.ip and (jl.2p are enhanced models in nonlinear acoustics, the Kuznetsov 

(1.5) Utt - bAut - c^Au = {^^{ut)'^ + |Vup)^ 
and the Westervelt equation 

(1-6) Utt - bAut - Au = {ji {l + ^) {ut)'^)^, 

are classical, well-accepted and widely used models governing sound propagation in fluids. As 
(11.11) and (11.21) . they are derived from the basic equations in fluid mechanics. The Kuznetsov 
equation is the more general one of these classical models, in particular the Westervelt equa¬ 
tion is obtained from the Kuznetsov equation by neglecting local nonlinear effects. Moreover, 
for a small ratio of z/ and Pr, that is, for small heat conductivity, dLSI) and (fTGl) can be 
regarded as simplifications of dni) and (11.21) . respectively. 

The classical models (|1.5p and (11.61) have recently been extensively investigated. In partic¬ 
ular, results on well-posedness for the Kuznetsov and the Westervelt equation with homoge¬ 
neous Dirichlet [KLn9| and inhomogeneous Dirichlet [KLV11| . [KL12| and Neumann [KLllj 
boundary conditions have recently been shown in an L 2 (ll)-setting on spatial domains 11 C M”' 
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of dimension n € {1,2,3}. Moreover, there are results on optimal regularity and long-time 
behavior of solutions for the Westervelt equation with homogeneous Dirichlet |MW11] and 
for the Kuznetsov equation with inhomogeneous Dirichlet |MW13] boundary conditions in 
Lp(D)-spaces where the spatial domain D C M” is of arbitrary dimension. 

On the contrary, mathematical research on higher order partial differential equations arising 
in nonlinear acoustics is still in an early stage. Well-posedness and exponential decay results 
for the homogeneous Dirichlet boundary value problems associated with (jl.ip and (na) in 
an L 2 (D)-setting where D C M”, n € {1,2,3}, have been shown in [Brul5] and [BK14] . 
respectively. In the present work we consider (jl.ip and (11.21) with inhomogeneous Dirichlet and 
Neumann boundary conditions in Lp(D)-spaces where the spatial domain D is of dimension 
n € N. We show global well-posedness and long-time behavior of solutions in an optimal 
functional analytic setting in the sense that the regularity of the solution is necessary and 
sufficient for the regularity of the initial and boundary data. While in |Brul5] and |BK14] 
the results were proved by means of appropriate energy estimates and the Banach fixed-point 
theorem, the techniques used in the present paper are based on maximal Lp-regularity for 
parabolic problems and the implicit function theorem in Banach spaces. 

We suppose that D C M”, n € N, is a bounded domain, i.e., an open, connected and 
bounded subset of the n-dimensional Euclidean space, with smooth boundary E. Let J = 
(0, r) for some finite T > 0 or J = M_|_ = (0, oo). We consider the inhomogeneous Dirichlet 
boundary value problem 


(1.7) 


(aA - dt){utt - bAut - c^Au) 
< (u, Au) 

{u,ut,utt) 


{k{utf + s\Vu\%t 

{9,h) 

{uo,Ui,U2) 


in J X D, 
on J X r, 
on {t = 0} X D, 


and the inhomogeneous Neumann boundary value problem 


( 1 . 8 ) 


(aA - dt){utt - bAut - c^Au) 
< {duU, dyAu) 

{u,ut,utt) 


{k{utf + s\Vu\%t 

{9,h) 

{uo,Ui,U2) 


in J X D, 
on J X r, 
on {t = 0} X D, 


where uq,ui,U 2 - D —>■ M and g,h : J xT ^ M. are given, u : J x D ^ M is the unknown, 
u{t,x), and a, 6, c and k are positive constants. Moreover, dyU = u ■ Va|r where u is the 
outer normal unit vector denotes the normal derivative of u. The parameter s € {0,1} allows 
us to switch between CH) and m- 

We point out that the present work extends the results from [BrulS] in several ways. 
First, while in |Brul5] the Blackstock-Crighton equation was considered with homogeneous 
Dirichlet boundary conditions, we also allow for inhomogeneous Dirichlet as well as Neumann 
boundary conditions. We are able to remove the restriction n € {1,2,3} on the dimension 
of the spatial domain D. Instead of L 2 (D), we consider (II.ip and (II.2p in Lp(D) where 
p £ (l,c)o) in case of the linearized equation and p > max{n/4 -|- 1/2, n/3} in case of the 
nonlinear equations (HZD and (jl.Sp . In particular, we require p £ (5/4, oo) in case n = 3 
and then p = 2 is admissible. Moreover, most notably, our conditions on the regularity of 
the data {g, h,uo,ui,U 2 ) are necessary and sufficient for the existence of a unique solution 
of the Blackstock-Crighton equation (within a certain regularity class/a certain subspace of 
Lp{J X D)). 

Our strategy for solving dLZl) and (|1.8p is to prove that their linearizations induce isomor¬ 
phisms between suitable Banach spaces and to apply the implicit function theorem. In some 
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sense these linearizations can be considered as a composition of a heat problem and another 
linearized problem for the Westervelt equation. While the linearized Westervelt equation can 
be handled similar as in [MWlll IMW13| . the heat equation has to be solved with higher 
regularity conditions. 

The paper is organized as follows. The purpose of Section[2]is to recall several facts we need 
on our way to global well-posedness and exponential stability of (11.71) and (II.Sp . In particular, 
we mention all function spaces we use, provide facts about the homogeneous Dirichlet and 
Neumann Laplace operator and list some important embeddings and traces. We also give a 
short review of the concept of maximal Lp-regularity for parabolic problems. Furthermore, 
we recall respectively prove optimal regularity results for the heat equation and the linearized 
Westervelt equation with inhomogeneous Dirichlet and Neumann boundary conditions. 

Section [3] is devoted to the inhomogeneous Dirichlet boundary value problem (|1.7p . First 
of all we consider the corresponding linear problem and represent it as an abstract parabolic 
evolution equation for which we show maximal Lp-regularity. This gives us optimal regularity 
for the homogeneous linear version of (II.7p . Based on the optimal regularity results for the 
heat and the linearized Westervelt equation from Section [2] we prove optimal regularity for the 
linear inhomogeneous Dirichlet boundary value problem in Proposition 13.51 The main result 
in this section is Theorem 13.61 which states global well-posedness of (jl.7p and immediately 
implies exponential stability (Theorem 13.7() . 

In Section 0] we treat the inhomogeneous Neumann boundary value problem (II.8|) . Here, 
we proceed analogously to Section [3l Theorem 14.81 provides local well-posedness for (jl.8l) . In 
the Neumann case, global well-posedness is shown for data having zero mean (Theorem 14.9p . 
Moreover, Theorem 14.111 states long-time behavior of solutions. 

In Appendix El we collect several facts about the homogeneous Neumann Laplacian. We 
outline how one finds a realization of the Laplacian with homogeneous Neumann boundary 
conditions such that its spectrum is contained in the positive half-line. 

In Appendix [B] we first study the temporal trace operator acting on a class of anisotropic 
Sobolev spaces. We present its mapping properties, provide a right-inverse and thus obtain 
its precise range space. Moreover, we construct functions with prescribed higher-order initial 
data. Second, we prove some so-called mixed derivative embeddings which are often used for 
checking the continuity of differential operators acting on anisotropic spaces. 

In Appendix Owe prove some higher regularity results for the heat equation with inhomo¬ 
geneous Dirichlet or Neumann boundary conditions and inhomogeneous initial conditions in 
a far more general framework than needed in the main text. In particular, we state explicitly 
all necessary compatibility conditions between initial and boundary data and show how they 
are used to contract a solution with high regularity. 

2. Preliminaries 

The purpose of this section is to introduce the notation we are going to use throughout the 
paper and to recall several important facts and results we need to prove global well-posedness 
and long-time behavior of solutions for (|1.7p and (jl.8|) . As already mentioned in Section [H we 
always assume that the spatial domain D C M”, re € N, is bounded and has smooth boundary 
T = dQ. We write J for a time interval and consider either J = (0, T) for some finite time 
horizon T > 0 or J = M_|_ = (0, oo). 

2.1. Function spaces, operators, embeddings and traces. The space BUC^{Q) con¬ 
tains all /c-times Frechet differentiable functions D —)■ M, whose derivatives up to order k are 
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bounded and uniformly continuous. For p G (l,oo), let Lp{Q) denote the space of (equiva¬ 
lence classes of) Lebesgue measurable p-integrable functions n —> M. We write W^{Q) for 
the Sobolev-Slobodeckij space and for the Bessel potential space of order m G [0, oo), 

where we have VF™(0) = if m G Nq. Moreover, W^{Q,;X) and H^{Q-,X) denote 

the vector valued versions and oWp{J]X) denotes the space of all functions u G Wp{J;X) 
with u(0) = 0. For p G [l,oo), q G [l,oo], s G M+, the Besov space Bp^y{Q) is dehned as 
(L„(n), n where m = \s] and (•, n indicates real interpolation. It holds that 

= W ;( n ) i!l £ R+ \N ah = ^5!) it p = , = 2. Moreovar, 

( 2 . 1 ) Bp g(fl) = (Wp{Q),W^{Q,))e^q, where 0 < k < s < m and s = (1 — Q)k + Qm. 

We always write X ^ Y the Banach space X is continuously embedded into the Banach 
space Y. Moreover, let L{X^Y) be the space of all bounded linear operators between X and 
Y. A linear operator A: X —>■ y is called an isomorphism if it is bounded and bijective. Then 
the closed graph theorem implies that A~^: Y —> X is also bounded and therefore A: X ^ Y 
is a homeomorphism. Now, let X and X be Banach spaces such that X Li^iodJ] X) where 
Li,ioc{J', X) is the space of locally integrable functions J —>■ X. For any a; G M we dehne the 
exponentially weighted space 

e'^X = {u€ Li^iociJ; X) : e-‘^^u G X}, 

equipped with the norm ||u||e“X = ||e““*u||x where e~‘^^u denotes the mapping [t i-A e~‘^^u{t)]. 

Let —Ad: T>{Ad) —?> Lp{Q), u e->■ —Art denote the negative Dirichlet Laplacian with 
domain 'D{Ad) = {u G Wp{Q): u = 0 on F}. We recall that the spectrum (t(—A^)) is 
a discrete subset of (0, oo) consisting only of eigenvalues = A„(—A^i), n G No with 
hnite multiplicity. We write \q > 0 for the smallest eigenvalue of —Ad- Moreover, the 
negative Neumann Laplacian —A^v : 'D{A]\f) Lp{Q), u i-a —Art with domain ©(A^r) = 
{u G Wp{Q): dyU = 0 on F} has a discrete spectrum (t(—A w) C [0, oo) which contains only 
eigenvalues = A„(—Aw), n G No, of finite multiplicity. Here, \q = 0 is an isolated point of 
(t(—A w) which can be removed when introducing the space Lpfi{^) = {u € Lp{Q): J^udx = 
0 } and considering — Aw,o: '^{Aj\f,o) —>■ Lp^^Xl) , u — Au with I?(Aw,o) = ^(^7v)FlLp^o(A^)- 
We then have (t(—A w,o) C (0, oo) where Xi = Ai(—Aw) > 0 is the smallest eigenvalue of 
—Aw,o- For details we refer to Appendix[Al 

We shall use the embeddings Wp{J) ^ BUC{J) and bFp(H) Wp{Q,) for s > t. We 
always write 7 _d = 'Ir and 'Jn = dy ■ \r = i' ■ (V-)|r for the Dirichlet and the Neumann trace, 
respectively. Moreover, 7 ^ = • |t=o denotes the temporal trace. Let B G {D,N}, Jd = 0, 
jw = 1- For p G (1, 00 ), /c G N and I G Nq the spatial trace 

w^+^{j-,Lpmnw^{j-,w^\n)) 

is bounded, see Appendix [Bj Furthermore, the trace 

(2.3) U^^BU: w;{n) ^ B;;J^-^/p(r) 

is bounded for every s G {Jb + 1/p, 00 ), cf. [Tri83[ Theorem 3.3.3]. The temporal trace 

(2.4) u\t=o : W^iJ- W;(L!)) n Lp(J; ^ 

is bounded for a G (1/p, 1], s G [0, 00 ) and s -|- 2 q: ^ N for a < 1. The same holds when the 
domain D is replaced by its boundary F. Finally, we mention that for p G (l,oo), G Nq, 
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r, cr G N, 0 € (0,1) we have the mixed derivative embedding 

(2.5) W^{n)) n W^{J; w^+^^in)) ^ J; 

where again ft can be replaced by F. For more general embeddings of this form we refer to 
Appendix |Bj 

2.2. Maximal Lp-regularity. Let J = (0,T) or J = = (0,oo) and assume p G (l,oo). 

We say that a closed linear operator A : T>{A) —>■ X with dense domain 'D{A) in a Banach 
space X admits maximal Lp-regularity on J if for each F G Lp{J;X) the abstract Cauchy 
problem 

(2.6) vt{t) + Av{t) = F{t), t e J, t'(O) = Vo, 

admits a unique solution u G ]E(J) = Wp{J;X) n Lp{J;T>{A)) for vq = 0. 

Furthermore, the abstract inhomogeneous Cauchy problem (12.61) is said to admit maximal 
Lp-regularity, if 

(2.7) (at + A, 7 i) :]E(J)^Lp(J;X) xtrE(J), v^{F,vo) 
is a homeomorphism. Then its inverse is the solution map 

(2.8) {dt + A,jt)~^ :Lp{J;X)xtr:E{J)^E{J), {F,vo)^v. 

If A : F’(A) —)■ X has maximal Lp-regularity on J, then the abstract Cauchy problem (12.61) 
has maximal Lp-regularity on J, cf. Section III.1.5 in [Ama95] . The following result is very 
useful and will be used several times throughout this paper. 

Lemma 2.1 (cf. |Ama95[ Proposition IILl.5.3]). Let a G M. Suppose that 

idt + a + A,^t)-W)^LpiJ-,X) xtrE(J) 
is a homeomorphism. Then 

{dt + A,-it)- e“E(J) ^e“Lp(J;A) x trE(J) 

is a homeomorphism. 


2.3. Optimal regularity results. In order to prove our results on optimal regularity for the 
linearized versions of (II.7p and (jl.8l) . we need optimal regularity results for the heat equation 
and the linearized Westervelt equation. We always let a, 6, c G (0,oo). 

2.3.1. Dirichlet boundary conditions. Recall that Xq > 0 always denotes the smallest eigen¬ 
value of the negative Dirichlet Laplacian in Lp(D). 


Lemma 2.2 f [LPS06[ Proposition 8]). Let p G (l,oo) and w G (0,aA^). 
boundary value problem for the heat equation 

' ut — aAu = f in M+ x D, 

(2.9) < u = 9 on M+ x F, 

u = uq on {t = 0} X D, 


Then the initial 


has a unique solution 

u G e-‘"ll„, EI„ = W;(M+; Lp{Ll)) n Lp(M+; W^{Ll)), 
if and only if the given data f, g and uq satisfy the regularity conditions 
(i) / G e-‘^Lp(M+ X D), 

(ii) no G 
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(iii) g G Mr = Lp{T)) n Lp{^+,Wp 

(iv) uo|r = g\t=o in the sense of traces. 


Lemma 2.3 f |MW13[ Lemma 5]). Suppose p G (l,oo), p 3/2, n G N and define ujq = 
min{6A^/2, /b}. Then for every u G (0, wq) there exists a unique solution 


W„ = W^{ 


.;Lp(12))nW3(M+;W^(f^)) 


•u G e Wu, 
of the linear initial boundary value problem 

utt — bAut — (?Au = /, in M+ x 0, 

(2.10) < n = g, on M+ x L, 

{u, ut) = (uo, ui) on {t = 0} X Q 


if and only if the data satisfy the following conditions: 

(i) / G e-‘^LpiR+ X n), 

(ii) uo G w^{Q), ui G 

(iii) g G e—Wr, Wr = Lp{T)) n W^iR+-W^-^^PiV)), 

(iv) g\t=o = tioir and if p > 3/2 also gt\t=o = i^i|r in the sense of traces. 


2.3.2. Neumann boundary conditions. Now we prove optimal regularity results for the heat 
equation and the linearized Westervelt equation with Neumann boundary conditions. Recall 
that Xi > 0 denotes the smallest eigenvalue of the negative homogeneous Neumann Laplacian 
in Lpfl(Q,). Let u = |n|“^ denote the mean of a function n : 12 —)> M and let g = 

|r|~^ '/r g dS for g : T ^ M. 

Lemma 2.4. Let p G (l,oo) \ {3} and oj G [0,oA)^). Then the inhomogeneous Neumann 
boundary value problem for the heat equation 

{ Ut — aAu = f in M_|_ x 12, 
duU = g on M+ X r, 
u = uq on {t = 0} X 12, 


admits a unique solution of the form u{t,x) = v{t,x) + w{t) with 

uGe ‘^M„^ 0 j (I^+i-1 ^p,o(1^)) I ~1 Lp(M+; PRp (12) n Lp^o(l^))) wt £ e ‘^Lp(M+), 

if and only if the data satisfy the following conditions: 

(i) /Ge-Lp(M+;Lp(12)), 

(ii) no G iyp“'/'’(12), 

(iii) g G e—M,, M, = iR+-, Lp(r)) n Lp(M+; Wp-^^^iT)), 

(iv) g\t=o = dvUo\r in the sense of traces if p > 3. 

If in addition f(t, •), no, g{t, ■) have mean value zero over 12 resp. L for all t, then w = 0. 

Proof. We first let to = 0. By [DHP031 Theorem 8.2], Lemma [4.51 and |Dor93l Theorem 2.4], 
the Nenmann Laplacian in Lpfi{n) with domain = T>[A]\j) n Lp^Qifil) has maximal 

regnlarity on R+. We therefore obtain a unique solution na G M„^o of the problem 

dtu^ — oAns = /a in M+ x 12, dyU^ = 0 on M+ x T, na(0) = 0 in 12 

for every given /a G Lp(M_|_; Lp^o(l^))- Furthermore, problem (12.lip admits at most one 
solution. Indeed, let us construct it as n = ni + n 2 + na where we first solve 


dtUi + /rni — oAni = 0 in M+ x 12, dyUi = g on M+ x T, ni(0) = no in 12, 
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for some sufficiently large /r > 0 with |DHP07l Theorem 2.1]. Next, we let U 2 solve the 
ordinary differential equation 

dtU 2 {t) = f{t) + U 2 { 0 ) = 0 , 

Finally, with /s = f — — we obtain U 3 as above. It is easy to check that v = U 1 +U 3 

and w = U 2 satisfy the assertion. The case cn > 0 can be reduced to the previous one by 
multiplying the functions u, f, g with 1 1 —>■ e“* and using that the spectrum of — oA^r^o + w is 
contained in ( 0 , 00 ). □ 

Lemma 2.5. Let p G (l,oo) \ {3} and u G (0,bA{^). Then the inhomogeneous Neumann 
boundary value problem 

{ utt — bAut = f in M_|_ x Ll, 

duU = g on R+ X T, 

{u,ut) = {uo,ui) on {t = 0} X n. 


admits a unique solution of the form u{t,x) = v{t,x) + w{t) with 

V € e-‘"W„,o, W„,o = w 2(M+; Lp,o(^^)) n W^{a) n Tp,o(f^)), 

W € Lpyoc{[0,oo)), Wtt € e““Lp(M+), 
if and only if the data satisfy the following conditions: 

(i) / G e-“Lp(M+ X n), 

(ii) no G Wp2(0), ni G 

(hi) g G e—w,, w, = Lp(r)) n W^{R+- ITp"^/'’(r)), 

(iv) g\t=o = 5jyno|r o,nd if p > 3 also gt\t=o = dyUi\Y in the sense of traces, 

/o°° = ^^“0 - ui ■ 

Moreover, w satisfies wttft) = f{t) + 6 |r||n|“^ 5 t(t) with w{0) = uq and wt{0) = hi. 

Proof. We start by proving sufficiency. First, note that dt : Wi, —> Hi/ is bounded and 

< 116^^*51^1/ which implies gt G e“‘^Hi/. Therefore, from 
Lemma 12.41 we obtain that the heat problem 

(ft — bAcp = / in x fl, = gt on M_|_ x F, </^(0) = ni in II, 

admits a unique solution of the form ip{t,x) = ipi{t,x) + P2{t) such that pi G e“‘^H„^o and 
dtip 2 £ e“‘^Lp(M+). In particular, since pi has zero mean over II, we have (/?2(0) = <^(0) = hi 
and <y9i(0) = ni — hi. For x G II and t G M+ we define n(l, x) = v{t, x) + w{t), where 

poo poo 

v{t,x) = — / (pi{s,x)ds and w{t) = — / ip 2 {s)ds. 

Jt Jt 

Clearly ut = ip, hence uu — bAut = / in H. Integrating the latter with respect to space, 
multiplying with |II|“^ and using the identity Audx = Jp d^udS, we deduce that w solves 
the ordinary differential equation 

wtt{t) = f{t) + 6|F||II|“^5t(t), w{0) = ho, wt{0) = hi. 

This implies vu — bAvt = / —/—6|F||II|“^5t in H. In what follows, we abbreviate v{t) = v{t, •), 
pi{t) = pi{t, ■) etc. and we let yK_(t) = 1 for t < 0 and XK_(t) = 0 for t > 0. Using to > 0, 
vt = Pi and the identity 

/ OO 
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together with Young’s inequality implies u G e Moreover, we have 

roo roo 

duv{t)\r =di,u{t)\r = - di,(p{s)\r ds = - gs{s) ds = g{t), 

Jt Jt 

roo POO 

vt{t) = ipi{t) = - dsipiis)ds = - {f{s) + bA(pi{s)-ds^p 2 is))ds 

Jt Jt 

f{s) ds + bAv{t) — ip 2 {t) 

and Vt{0) = ¥^ 2 ( 0 ) = ui — tti. Altogether, 6Au(0) = Ut(0) + rti + f{s) ds = bA{uQ — uq) 
and 5jyu(0)|r = g{0) = du{uo — uq) which implies that u(0) = uq — uq in Q. 

To verify necessity of (i)-(v), we assume that u{t, x) = u(t, x) + w{t) with v G e“‘^W„^o and 
wtt G e“‘^Lp(M+) is a solution of ()2.12p . We have f = e^^'Vu + + bA{e^^vt). Since 

e^^ut = {e^^v)t — ioe‘^^v G Lp(M+; Wp{Q) D Lpfl{Q,)) and e‘^^Vtt = — 2u}e‘^^Vt — w^e^^u G 

Lp(M+; Lp^o(f^)) we conclude e‘^^f G Tp(M+ x Q) and (i) is verified. Concerning (ii) note 
that exponential weights do not affect the initial regularity. Due to hTp(M+) i?f7C'(M+) 
we infer v\t=o G Wp{Q) D Lp^o(f^)) hence u|t=o = u|t=o + w^|t=o £ Furthermore, 

we have G Wp (M+; Lp^o(^)) C Lp{M^;Wp{i}) n Lp^o(^))- Applying the temporal trace 

()2.4p with a = \ and s = 0 implies vt\t=o G Bp^p^ = Wp hence ut\t=o = 

vt\t=o + wt\t=o G VFp In order to check (iii), we apply the spatial trace (|2.2p with 

A: = / = 1 to e‘^*u G Wu,o which gives us e‘^^di,u\r G as claimed. Next, note that 

(12.,ip applied to uq G Wp{Q) implies di^UQ G i?p,p^'^^(r) = VFp ^^^(F). Moreover, we have 
g G IFpi(M+; IFp"^/^(r)) -A BUC{^+, IFp"^/^(F)), hence S^uolr = g\t=o in ITp“^/^(r). From 
ui G Wp we obtain d,yUi G i?p,p^'^^(F) = IFp ^^^(F) for p > 3 and, from gt G Bu, 

using p2.4p with a = 1/2 — l/2p and s = 0, we get gt\t=o £ Bp^^^^iT) = Wp ^^^(F) if 
p > 3. Altogether, dyUi\Y = gt\t=o in Wp for p > 3. For (v), note that integrating 

utt{t) — bAut{t) = f{t) with respect to time yields ut{t) — bAu{t) = — f{s) ds, hence in 
particular ut(0) — bAu{0) = — f{s) ds. Therewith the proof of necessity is complete. 

Finally, we show that (I2.12p has at most one solution. To this end, suppose we have given 
two solutions of p2.12p . Their difference u solves 

utt — bAut = 0 in M+ x D, dyU = 0 on M_|_ x F, u(0) = ut(0) = 0 in 11. 
Furthermore, u = ut solves the heat problem 



lit — bAuf = 0 in M+ x 11, di,u = 0 on M_|_ x F, u = 0 in D, 

which implies ut = 0. Hence u is constant which together with {t(0) = 0 implies u = 0. □ 

Lemma 2.6. Let p G (l,oo) \ {3} and cjq = min{6A^/2, c^/6}. Then for every oj G (0,wo) 
the initial boundary value problem 


(2.13) 


Utt - bAut - c^Au = /, 

< dyU = 0 , 

{u,Ut) = (uo,Ul) 


in M+ X H, 
on M+ X F, 
on {t = 0} X H, 


has a unique solution 

u G e“‘^W„,o, 


,,0 = w^{R+, Lp,o(ii)) n w/(M+; w^{n) n Lp,o(ll)) 
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if and only if 

(i) /ee-Lp(M+;Lp,o(^^)), 

(ii) uq G Wp{il) nLp_o(^) 0 'iT‘d ui € Wp nLp^o(f^) such that d^uo\r = 0 and if p > 3 

additionally 9jyixi|r = 0. 

Proof. In [MWll] the result was established with Dirichlet instead of Neumann boundary 
conditions. Here, we just point out the main steps of the proof. We represent (I2.13P as 



Let us consider the space X = P(A7v,o) x -bp,o(^) and the operator A: V^A) X given by 

(2.14) = 


First we show that there is some n > 0 such that v + A admits maximal Lp-regularity on M_|_ 
by a perturbation argument. Choosing a decomposition of ^4, ^4 = ^1 + 742 with 


(al -I \ 
V 0 —bXjsifl J 


and A 2 


( —al 0 \ 

I^-c^Aat^o oJ 


for some a > 0, it turns out that operator Ai: 11{A) X admits maximal regularity on M_|_ 
due to Lemma [23] and the fact that the bounded operator {dt + a) \ oVFp (M+; ^(AAr^o)) 
Lp(M+;D(AAr,o)) is invertible. Moreover, since A^'. X ^ X \s bounded, on the strength of 
Proposition 4.3 and Theorem 4.4 in |DHP03| there exists some z/ > 0 such that v -\- Ai + A 2 
admits maximal regularity on M+. 

Next we claim that the spectral bound s(—A) = sup{Re(A): A G a{—A)} of —A is given by 
s(—= —ujQ. This follows analogously to [MWlll Lemma 2.4] if one replaces the Dirichlet 
eigenvalues by the ones of 

Since s(—^4) = —cjq < 0, the spectral bound of —A + uj equals uj — ujq which is strictly 
negative as long as a; G [0,tLio). Hence, for each uj G [0,Ci;o) the operator A — uj has maximal 
Lp-regularity on M+ by |Dor93l Theorem 2.4], that is 

(dt+A- UJ, Jt): W^pHK+; n Lp(R+;V(A)) ^ Lp(M+; A) x (A, D(i))i_i/p,p 
is an isomorphism. Employing Lemma (2. II we conclude that for every uj G [OjCJq) the operator 
(dt + A,jt): e-‘^(WpHu+;X) n Lp(R+;V(A))) ^ e-‘^Lp(R+;X) x (X, D(A)),_,/p^p 
is an isomorphism. 

It is easy to check that {u,ut) G e“^(ITp (R+; A) n Lp(M+;D(A))) implies u G e~‘^Wu,o- 
Moreover, we have / G Lp(M+; Lpfi{Q)) and uq G P(A 7 v), i.e. uq G ITp (D) nLp^o(^) such that 

duUolr = 0. Finally, ui G {Lp^o{n),'D{AN))i-i/p^p = Wp“^^^(D) n Lp,o(^^) with d^ui\r = 0 
where the trace exists if p > 3. This concludes the proof. □ 


Finally we arrive at our optimal regularity result for the linearized Westervelt equation 
with inhomogeneous Neumann boundary conditions. 
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Lemma 2.7. Letp^ (l,oo)\{3} and set uq = mm{bX^/2,c^/b}. Then for every uj ^ (0,wo) 
the linear initial boundary value problem 

utt — bAut — (?Au = /, in x 

(2.15) < duU = g, on M+ x F, 

{u,ut) = {uo,ui) on {t = 0} X Q, 
admits a unique solution of the form u{t,x) = v{t,x) + w{t), where 


V G e-“W„,o, W„,o = Wp2(M+; Lp,o(f^)) n W^{n) n Lp,o(^^)), Wu G e--Lp(M+) 

if and only if the data satisfy the following conditions: 

(i) / G e--Lp(M+ X n), 

(ii) 1X0 G w^iQ), ui G 

(iii) g G W, = Lp(r)) n W^{R+-Wp-^^^iT)), 

(iv) g\t=o = 5jyUo|r o,nd if p > 3 additionally gt\t=o = dyUi\Y in the sense of traces. 


Proof. From Lemma 12.61 we obtain uniqueness. In order to show necessity of (i)-(iv) one 
proceeds as in the proof of Lemma 12.51 It therefore remains to show sufficiency. Let 5 > oj. 
From Lemma 12.51 we obtain that 

{ ptt - hAipt = f - fs in M+ X 12, 
d^(p = g, on M+ x F, 

{(p, p>t) = (uo, Ui) on {t = 0} X n, 

where fs = f{s) ds + ui — bAuo), admits a unique solution (p{x, t) = ipvix, t) + piw{t) 

such that p>y G e^'^W^^o and G e“‘^Lp(M+). Next, Lemma YIM implies that 

- hAOy^t - = /^ - /^ + - c^Ap., in R+ x 12, 

< d,y6y = 0, on M+ X F, 

{0v,0v,t) = (0,0) on {t = 0} X 12, 

has a unique solution 6 G e“^Wn,o- Furthermore, we define 9w as the solution of the ordinary 
differential equation 

9wpt{t) = c^A(py{t) + fs{t), 6'u;(0) = 0, = 0. 

Then v = ipy + 6.,] and w = pw + Oy, satisfy the assertion and we are done. □ 


Remark 2.8. If we consider (12.151) on a finite time interval J = (0, T) instead of M_|_, we 
may set a; = 0 and obtain a unique solution u G Wp{J; Lp{Q,)) n Wp{J] VF^ (12)) if and only if 
conditions (i)-(iv) (with u; = 0) hold. 


2.4. Analysis in Banach spaces. For later use in the proof of global well-posedness of 
(11.71) and (|1.8p we will now recall the concept of analytic mappings in Banach spaces and the 
analytic version of the implicit function theorem. The remainder of this section is collected 
from Section 15.1 in [Dei85j . 

Let X and Y be Banach spaces over the same field IK = M or K = C and let 1/ C X be 
open. Then F: U ^ Y is called analytic at xq G C/ if there is some r > 0 and continuous 
symmetric fc-linear operators Fj^: = X x ■ ■ ■ x X ^ Y for k > 1 such that 


k=l 


and F{xo + h) = F{xo) + ^ Fk{h^). 

k=l 


< oo 
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for h & X, ||/i|| < r. Here, ||Ffc|| = sup{||Ffc(/i^)||: ||/i|| < 1}. We then necessarily have 
Fk = The map F: U is called analytic if F is analytic at every xq G U. In 

particular, every bounded linear map F : X ^ Y is analytic. 

Theorem 2.9 (Implicit Function Theorem, cf. [Dei85[ 15.1]). Let X, Y and Z be Banach 
spaces over the same field IK = M and IK = C. Assume U C X and V CY are neighborhoods 
of Xq £ X and uq G Y, respectively. Furthermore, suppose 

(i) F: U xV^Z, {x,y) i-A F{x,y) is continuous, 

(ii) the Frechet derivative F^: U xV ^ L{X,Z) of F with respect to x is continuous, 

(hi) F{xQ,yo) = 0 and Fx{xQ,yQ)-. X ^ Z is an isomorphism. 

Then there exist balls B^{xq) C U and Hp(yo) C V and a unique map (p: Hp(?/o) —>■ Bn{xQ) 
such that ip{yo) = xq and F{ip{y),y) = 0 for all y G Hp(yo). The map ip is continuous. If 
furthermore, F is analytic, then ip is analytic in some neighborhood of yQ, in particular on 
some (possibly smaller) ball B^iyo) C Bp^yo). 


3. The Dirichlet boundary value problem 

In this section we prove global well-posedness and exponential stability for ()1.7p . First of all, 
we consider the linearized version of the inhomogeneous Dirichlet boundary value problem and 
represent it as an abstract evolution equation. We show that this abstract equation admits 
maximal Lp-regularity and derive an optimal regularity result for the linearized equation 
associated with (HZl). Then we use the implicit function theorem to construct a solution of 
the nonlinear problem ()1.7p . Exponential decay of this solution is an immediate consequence. 


3.1. Maximal Lp-regularity for the linearized equation. Suppose J = (0, T) or J = R+. 

For / G Lp{J X D) we consider the initial boundary value problem 

(oA — dt){utt — bAut — (?Au) = f in J x D, 

(3.1) < {u,Au) = {g,h) on J X F, 

{u,ut,utt) = iuQ,ui,U 2 ) on {t = 0} X D, 


where uq,ui,U 2 '. D —M and g,h'. J x T —>■ M are the given initial and boundary data, 
respectively. In order to address the problem of maximal Lp-regularity for the linearized 
equation, we represent (13.ip with g( = /i = 0 as an abstract Cauchy problem 


\ \ 0 0 -anAj J \utt - bAout - c^ Aud ) \-fJ 

This motivates us to consider the Banach space 

(3.2) X^ = VHAd?) X V{Ad) X Lp(D) 

and the densely dehned linear operator : T>{A^) —>■ X^ given by 


/ 0 -I 0 \ 

(3.3) A^ = -c^Ad -bAjo -I , V{A^) = P((Ab)") x VHAnf) x V{Ad). 

y 0 0 —aAjo j 

Therewith, we may write (|3.ip as an abstract evolution equation 

dtv^ + Av^ = F, u^(0) = u^ 
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if we define 
(3.4) 


v^ = 


u 

Ut 


- 

Un — 


UO 
Ui 

.U2 — h/S.£)Ui — , 



^utt - bAout - c^AdUj 
First of all we will treat the issue of maximal Lp-regularity of : T>{A^) —>■ on M+. 


Proposition 3.1. Let p € (1, oo). There is a constant p > 0 such that p + A^ has maximal 
Lp-regularity on M_|_. 

Proof. Let a > 0. We decompose A^, A^ = Af + A 2 , where 

(al -I 0 \ / -al 0 0\ 

Ai = j 0 —bAo —I j and A 2 = I —c^An 0 0 j . 

\0 0 -aAol \ 0 0 0 / 

First we show that A^ : 'D{A^) —> has maximal Lp-regularity. To this end, we consider 

the Cauchy problem vt + Afv = F, vo = 0 and show that for each F G Lp(M+;X^) there 
exists a unique solution v € 1F/(M+;X^) n Lp(R+;P(^^)). With v = {vi,V 2 ,V 3 )~^ and 
F = (/i, / 2 , fsV, we explicitly have 

dtvi + avi -V 2 = fi, ui(0) = 0, 

dtV 2 - bADV 2 - V 3 = f 2 , U2(0)=0, 

dtvs - aAoVs = /s, U3(0) = 0. 

Since we know from Lemma 12.21 that the homogeneous heat equation admits maximal Lp- 
regularity, we obtain that for all /a G Lp(M+ x Ll) there exists a unique solution 

U3 G W/(M+;Lp(0))nLp(M+;P(Az,)). 

Moreover, as /2 + U 3 G Lp(M+; ©(A^))), Lemma IC.51 implies that there is a unique solution 

V 2 G VF/(M+;P(Az 5 )) nLp(M+;P((Az 5 ) 2 )) 

Now, note that for a > 0 the operator (dt -t- a ): o 1F/(R+; L>((A£))^)) Lp(M+; P((A£))^)) is 
invertible. Since fi + V 2 ^ Lp(M+; P((A£))^)) we obtain a unique solution 

vi{t)=[ 112 ( 5 )) ds, 

Jo 

which satisfies vi G PF/(M+;^((Ai:))^)). Altogether, we conclude that Ai : V{A^) —)> X^ 
admits maximal Lp-regularity. 

Moreover, by the fact that A 2 ■ X^ X^ is a bounded linear operator. Proposition 4.3 
and Theorem 4.4 in [DHPr)3] imply that there exists some p > 0 such that p + Af -|- has 
maximal regularity which concludes the proof. □ 


In order to show maximal regularity for the operator A^ : T>[A^) —>■ X^, we need the 
following result on its spectrum. 

Lemma 3.2 (cf. |BK14[ Lemma 3.10]). The spectral bound s{—A^) = sup{Re(A): A G 
a{—A^)} of —A^ is given by s{—A^) = —ujq, where Uq = min{aA^, 6A^/2, c^/ 6 }. In 
particular, i/Re(A) < ojq , then A G p{A^). 
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Theorem 3.3. Let p € (l,oo) and u € [0,a;o) where oJq = min{aA^, 6A^/2, c^/6}. Then 
: T>{A^) —>• has maximal Lp-regularity on M_|_ in the sense that 

{dt + A^,^t)--e-‘^{W^{R+,X^)nLp{R+,V{A^))) 

^e-‘"Lp(M+;X^) x {X^,V{A^)),_yp^p 

is an isomorphism. 

Proof. We follow the proof of Theorem 2.5 in |MW11| . From Proposition 13.II we know p + A^ 
admits maximal regularity on M+ for some p > 0. Multiplying + A^v^ = F hy e“^* shows 
that A^ has maximal Lp-regularity on bounded intervals J = (0,r). Lemma [3.2l tells us that 
spectral bound s{—A^) = —ujq of —A^ is strictly negative. Hence s{—A^ + oj) = u — ujq < 0 
as long as w G [ 0 ,a;^). From [Dor93l Theorem 2.4] we deduce that A^ — uj admits maximal 
Lp-regularity on M_|_ for every a; G [ 0 , 0 ;,^), that is, 

{dtPA^-uj.yt)-. W;(M+;X^)nLp(M+;P(H^)) ^Lp(M+;W^) x P(H^))i_i/p,p 

is an isomorphism. Now Lemma l2.II implies the result. □ 


Corollary 3.4. Let p G (l,oo) \ {3/2} and define ojq = min{oA^,6A^/2,c^/6}. Then for 
every oj G (OjCJo) the linear Dirichlet boundary value problem 


(3.5) 


{aA — dt){utt — bAut — (?Au) = f in M+ x H, 

< {u, Au) = (0,0) on M+ x F, 

{u, ut,uu) = {uq,ui,U 2 ) on {t = 0} X H, 


admits maximal Lp-regularity in the sense that there exists a unique solution 
u G e-‘"E„, E, = W3(]R+; Lp(H)) n W^{R+; W^{n)), 

if and only if 

(i) / G e--Lp(M+ X n), 

(ii) uo G VFp(H), ui G kFp U2 G VFp with uojr = = ^^i|r = 0 and, if 

p > ‘ij'l, also Auijr = U2|r = 0 in the sense of traces. 


Proof. Based on the choices of X^ and 'D{A^) in ()3.2p and (13.31) . it is straightforward to 
check that the condition G e“‘^(kFp (M+; D Lp(M+; L>(H^))) where is given by 
()3.4p implies u G e“^(Eu n kFp (M+; kFp (H))). Since the mixed derivative embedding gives 
us Eu ^ VFp (1R+; Wp (H)), we arrive at u G e^^E^. Next, we determine {X,V{A))i_iip p. 
It is trivial that (L>((A£))^),P((A£))^))i_i/pp = P((A£))^), i. e. we have uq G Wp{Q.) with 
u|r = Attjr = 0. Moreover, since for p G (1, oo) we have 2/p G M \ N unless p = 2, ()2.ip gives 
us 


{w^m,w^{n)fi_yp^p = B^p-^/p{n) = Wp"-2/p(H), 
(Lp(L!),iyp2(L!))i_i/p,p = B^pf/P{n) = Wp2-2/P(H). 

Moreover, interpolation with boundary conditions as in [AmaOQl Section 4.9] yields uijr = 
Auijr = U2 — bAjjui — c^A/juojr = 0. Hence, we have uq G PFp(H), ui G Wp and 

U2 — c^Aduq — bAoUi G kFp which is equivalent to uq G Wp{Ll), ui G kFp and 

U2 G VFp The result now follows from Theorem 13.31 □ 
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We now arrive at the final result for this section and prove optimal regularity for the linear 
initial boundary value problem (j.S.lIl in the sense that the regularity of the data /, g, h, uq, 
ui and U 2 are necessary and sufficient for the existence of a unique solution u G 

Proposition 3. 5. Let p £ {l,oo)\ {3/2} and define coq’= mm{aXQ ,bXQ /2, cP/b}. Then for 
every oj £ (0,a;o) the linear initial boundary value problem 


(3.6) 


(aA - dt){utt - bAut - c^Au) 
{u,ut,utt) 


f in M_|_ X n, 

{g, h) on R+ x F, 

(uO) '*^1)'*^ 2 ) on {t = 0} X n, 


has a unique solution 

u £ e-‘"E,, = W3(M+; Lp(0)) n W^{R+- W^{n)), 


if and only if the data satisfy the conditions 
(i) / G e-“Lp(M+ X n), 

(ii) uo £ W^{n), ui £ U2 £ 

(hi) g £ e-Eg,r, Eg,r = Lp{V)) n W;(M+; Wp""^/^(r)), 

h £ e—Lp(r)) n W^iR+-W^-^/P{T)), 

(iv) uoir = 514=0, -wilr = 5t|4=o, Auo|r = h\t=o and, if p > 3/2, also Ani|r = ht\t=o and 
^2|r = 54t|4=o hold in the sense of traces. 

Moreover, the solution fulfills the estimate 


ll“lle-“E„ ^ ll/lle-“Lp + ll5lle-“Fg,r + ll^lle-“F,i,r + ll'“o||w4 + ||R1 ||^^4-2/p + ||«2 ||^2-2/p. 

Proof. First we show necessity of (i)-(iv) for the existence of a unique solution u £ e^^^E^ 
of (j3.6p . In the proof of Corollary 13.41 we already mentioned that E^ iyp(M+ x 0). 
Since {e^^u)u = oj‘^eP^u + 2u]e‘^^ut + eP^uu and {e^^u)tu = + 

3ijj‘^e^^ut + 3uje^^uu+cP^utu-i the assumption that u £ e^^^E^ implies eP^ut £ Lp(M+; VFp (n))n 
W;(M+; W2(fl)) n W^{R+-Lp{Ll)), eP^uu £ Lp{R+-W^{Ll)) n W;(M+; Lp(fl)) and e^^utu G 
Lp(M+ X Ll), hence e^*/ = e^*(aA — dt){uu — bAut — (fiAu) £ Fp(M+ x Lt) and (i) follows. 
Next, we show (ii). The embedding VFp(J) BUC{J) implies uq G Wp{Ll) whereas the 
temporal trace (12.41) with a = 1, s = 2 and a = 1, s = 0 gives us the desired regularities of 
ui and U 2 , respectively. 

For u £ e““Et( the spatial trace (|2.2p with k = 1 and I = 2 implies ir|r = 5 G e“‘^Fp^r and for 
Au £ e“^W„ the choice k = I = 1 gives us Aiijr = h £ e^^^F/i^r- This shows (iii). 

Using VFp(J) “—)> BUC{J), the spatial trace (|2.3p . the temporal trace (|2.4p and the mixed 
derivative embedding (|2.5I) one shows (iv). We have 

uoir = 5|t=o in W^-^Ip{V), AuqIt = h\t=^ in W^-^/p{V), 

Rilr = 5t|4=o in B^~^/p{T), Aui|r = ht\t=Q in BI~^/p{T) if p > 3/2, 

R2|r = gtt\t=o in bI~^/p{T) \i p > 3/2. 

It remains to show that conditions (i)-(iv) imply the existence of a unique solution u £ 
e^^^E^ of (|3.5I) . Since we are dealing with a linear partial differential equation with constant 
coefficients, we may interchange the order of differentiation on the left-hand side and consider 
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the subproblems 


(3.7) 


f 2 

wtt — bAwt — c Aw = f 
< w = ah — gt 

{w,wt) = {aAuQ — ui,aAui — U 2 ) 


in M_|_ X n, 
on M_|_ X r, 
on {t = 0} X n, 


and 

(3.8) 


' aAu — ut = w in M_|_ x 0, 

< u = 5 on M+ X r, 

u = uq on {t = 0} X 0. 


From condition (ii) we obtain aAuQ — ui G Wp{Q) and aAui — U 2 G Wp ^^^(fi). Furthermore, 
(iii) implies ah — gt G e“‘^Wr. On the strength of Lemma [2l3] we obtain that (|3.7p admits 
a unique solution w G e~‘^Wu- Now we use Corollary 1C.41 with I = 1 and A: = 2 to solve 
(|3.8p and obtain that (|3.8p has a solution u G This concludes the proof of sufficiency. 

Uniqueness follows from Corollary 13.41 □ 


3.2. Global well-posedness and exponential stability. Based on ProDosition l3.5l we now 
show that there exists a unique global solution of the nonlinear initial boundary value problem 
()1.7p which depends continuously (in fact, even analytically) on the (sufficiently small) initial 
and boundary data. Moreover, we prove that the equilibrium u = 0 is exponentially stable. 

Theorem 3.6 (Clobal well-posedness - the Dirichlet case). Let p > max{n/4-|- 1/2, n/3}, 
p 7^ 3/2 and define Uq = min{aA^, /2, c^/6}. Suppose 

uo G Ui G U2 G 

(3.9) g G e-‘"F3,r, ¥g,r = Lp{T)) n 

h G e-‘"F;,,r, = W^-^/^p{R+-L p{V)) n W^{R+-, lUp"-'/^(F)). 

with Uoir = g\t= 0 , «i|r = gt\t=o, Auo|r = h\t=o and, if p > 3/2, also Aui|r = ht\t=o, 
U 2 \v = gtt\t=Q- 

Then for every oj G (0,a;o) there exists some p > 0 such that if 

llfi'lle-“Fg,r + ll^lle-“Fh,r + ll''^o||w/ + 1^1 lly[^4-2/p + ||u2||^2-2/p < p, 
the nonlinear initial boundary value problem (HZD admits a unique solution 

(3.10) u G e-‘"E„, E„ = n bN(]R+; W^{Ll)) 

which depends analytically on the data (ITOI) with respect to the corresponding topologies. 
Moreover, conditions ()3.9p are necessary for the regularity of the solution given in ()3.10p . 

Proof. Employing the results on the linearized problem (|3.ip from Section 13.11 we will now 
construct a solution of the nonlinear initial boundary value problem (|1.7p which we linearize 
at u = 0. Hence, the solution will be of the form u = u* + tt,, where u* solves the linearized 
problem (|3.ip for the data (/ = 0, g, h, uo,ui,U 2 ) and u, satisfies homogeneous boundary and 
initial conditions. We will hnd the (small) deviation u, from u* by application of the implicit 
function theorem to the map 

G: e~‘^Ru,h x e“‘^E„ e~‘^Lp{R+ x H), 

(u,,u*) i-A D{dt,A)u, - {k{{u, + u*)t)^ - s\V{u, + u*)P)u 


(3.11) 
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where the differential expression D{dt, A) is given by D{dt, A) = (aA — dt){df — bAdt — c^A) 
and = {u £ Ku'. u{0) = ut{0) = utt{0) = 0, ix|r = An|r = 0}. Explicitly, we have 
D{dt, A)u, = + (a + b)Au,^tt + <?Au,^t - abA^u,^t - ac^A'^u,. 

Step 1: The implicit function theorem applies. First of all, we will now verify the assumptions 
of the implicit function theorem (Theorem 12.91) . 

Step 1(a): G is analytic. The mixed derivative embedding (jB.lah implies that the linear 
maps 


u, i-A ■ 

u, i-> Au,^tt : 
u, !->■ Au,^t '■ 
u, i-A A^u,^t'- 
u, i-A A^u ,: 


e ^ e ^Lp(M_|_ x n), 

e e ^Lp(M_|_ x n), 

^ e-‘"(Lp(M+; W^{^)) n iy;(M+; Lp{^))) 
e e ^Lp(M_|_ x Q), 

e-‘^Eu,h^e-‘^W^{R+-,Lp{n)), 


are bounded and therefore analytic. Hence n, D{dt, A)u,: e e ‘^Lp(M+ x H) is 

analytic. Next, note that for p > 1/2 + n/4 the embedding 


E„ -A it;(m+; w^in)) n w2(m+ x n) 

-A H^+^/P+%R +--A BUC\R+]BUC{n)) 


holds. In particular, it holds if e > 0 is sufficiently small and 4 — 2/p — 2e — n/p > 0. Such an 
e > 0 exists if p > (n + 2)/4. Moreover, on the strength of the mixed derivative embedding 
theorem and the Sobolev embedding theorem we conclude similar as in the proof of Lemma 


6 in [MW13| that 


E„ -A ITp2(M+; w^in)) n IT3(M+; Lp{n)) 

-AlT|p(M+;W-2-2e(L!)) 

-AlT|p(M+;L2p(f^)) 


for 0 — I € [0,1] and e > 0, 

for e > 0, 

for 0 > ^ + e, 

for 0 < 1 - 


provided e > 0 is sufficiently small and p > 1/2 + n/4. Furthermore, we observe that 
e“2^Lp(M+ X H) e““Lp(M+ x H) since e^* < e^^* for a; > 0. Prepared like that, we 
estimate 


Wftgthp < II /til Lap II /t II Lzp < ll/llE„||ff||E„, 

\\{ftgt)t\\Lp < ll/n||Tpll5i||Loo + ll/dlToollfftt||Tp < II/IIe„||5||e„, 

II (/t<?t)n||Tp A ll/nt IIUp 11^4 11 Too T ^ll/ti IIt 2 p \\gtt IIt 2 p ff II/til Too ll*?in||Tp ^ II/IIe„ IlffllEu ) 

and conclude that {f,g) i-A ftgt- E^ x E„ ^ IFp (M+; Lp(H)) is bilinear and bounded, thus 
analytic. Setting u; = n, + n* in 

(3.12) 

and choosing / = e'^^n, and g = e^^^n* proves that 

(3.13) (li., tt*) (((tt, + ■ e~‘^'^u,h X e““E„ ^ e““Lp(M+ x H) 
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is analytic. It remains to show analyticity of the map 

(rt,, u*) I—)■ (I V(rt, + tt*): e x e —>■ e ‘^Lp(M_|_ x 11). 

Note that we have the embeddings 

E« W^{n)) -A BUC{R+-, W^{n)) ^ BUC{R+] BUC^{n)), p > n/3 

E« iyp^(M+; Wp^(II)) n %^(M+ X 0) -A H^+^/^P+^R+-,H^-^/P-‘^^{n)) 

-A ii^ 2 V(e+; H^-^/p-^^{n)) -A ii^ 2 V(e+ x ii), p> n/e +1/3 

Therewith, we obtain the estimates 

||V/-V^IIl, < \\s/f\\L,Jgh2, < ll/lkll^lk, 

||(V/ • Vg)t\\L, < ||(V/)t||L,J|V5||L,, + ||V/|U,J|(V5)t||L,, < ||/lkll5llE„, 

||(V/ • V<7)it||L, < ||(V/)«|UJ|V5 ||l^ +2||(V/)t||i,J|(V5)i||L,, + ||V/||i^||(V<7)it||L, 

< II/IIe„||5'IIe„, 

and conclude that (/, g) i-A V/ • S/g : E^j x E^ —>■ ITp (M+; Lp{Q)) is bilinear and bounded, thus 
analytic. Moreover, we have 

= ((e‘"ku;)2)ii -4a;((e‘"ku;)2)i + 4a;2(e“ku;)2. 

By setting tc = u, + tt*, / = e'^^u, and g = e‘^^u- 1 , we are done. Altogether, we have that 
G: e~‘^Eu,h x e^^^E^ —>■ e“‘^Lp(M+ x fl) is analytic. 

Step 1(b): Du,G{D,Q): e“^Lp(]R+ x fl) —>■ e~^Eu,h is an isomorphism. The Frechet derivative 
of G with respect to u, at (0,0) is given by 

Du.G{0, 0)[u] = (uA - dt){utt - (?^u - bAut). 

The map Du,G{0,0): e“‘^Lp(M+ x fl) —>■ e~‘^Eu^h is an isomorphism since, according to 
Corollary 13.41 for every / G e~‘^Lp{R^ x II) the equation (aA — dt){utt — c^Au — bAut) = / 
admits a unique solution u € e~‘^Eu,h- 

Step 2: Construction of the solution. On the strength of the Implicit Function Theorem 
there exists a ball Bp(0) C e^^E^j with sufficiently small radius p > 0 and an analytic map 
p: Bp{0) C e^^^Eu —>■ e~‘^Eu,h, u* i-A tt, = p(u*) satisfying <p(0) = 0 and G(<p(ti*), u*) = 0 for 
all tt* G Bp{0). Hence, whenever u* satisfies the boundary conditions u*|r = g, Au*|r = h 
and initial conditions u*|t=o = uq, Ui,^t\t=o = '“i, Ui,^tt\t=o = ^2 which is the case if we 
define u* G e^'^E^j to be the unique solution of (I3T]) with (/ = 0 ,uo,ui,U 2 , g, h), then 
u, + u* = <p(u*) + u* solves (II.7p . 

Step 3: Dependence of the solution on the data. It remains to show that the solution u G 
e^^^Eu depends analytically on (p,/i, uq, ui, U2). To this end, we dehne the spaces 

1 := e-^^Ep X e-^'E,, x W^{D) x W^-‘^/p{D) x W^-^/p{D), 

D := {{g,h,UQ,ui,U 2 ) G D: uo|r = 5|t=o,«i|r = 5t|i=o,'«2|r = 51*11=0 if P > 3/2, 

Auoir = /i|i=o, Aitilr = ^i|i=o if P > 3/2}. 

From Proposition 13.51 with / = 0 we obtain that u* depends linearly and continuously and 
thus analytically on (p,/i, uq, ui, 1x2) G B. Moreover, u* i-A tt, = p(ti*) is analytic on Bp(0) 
and therefore tt, G e~‘^Eu,h depends analytically on the data (p, h, uq, ui, U2) & B. Altogether, 
u = ti, + u* enjoys the same property which concludes the proof. □ 
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An immediate consequence of Theorem 13.61 is that the global solution u G e of ()1.7p 
decays to zero at an exponential rate. 

Theorem 3.7 (Exponential stability - the Dirichlet case). Under the same assumptions as 
in Theorem ESI the solution u decays exponentially fast to zero as t ^ oo, in the sense that 

\Wii)\\wf + \Wtit)\\^i- 2 /p + \\utt\\^, 2 - 2 /p < Ce t>0, 

for some C > 0 depending on the boundary and initial data g, h, uq, ui and U2- 

Proof We have u G e-"^WpH]R+; W^{n)) -A e-^BUC{R+; W^{n)), hence 

u G BUC{R^,Wp{Q)), ^Cie with Ci = ||e‘^ u\\Buc{K+-,wf)- 

Furthermore, G ]HI(M+) ^ Wp for j G {0,1,2}. Therefore, we obtain 

ut G BUC(R+, ||ut(t)||^ 4 - 2 /p < C 2 e-“* with C 2 = lls^c(R+;lUp"-"/P)- 

Finally, from utt G ]H[(M+) we deduce that 


utt G BUC{R+, W2-2/P(0)), \\uttit)\\^ 2 - 2 /p < C 3 e-‘"‘ with C 3 

and the claim follows. 


□ 


4. The Neumann boundary value problem 

In this section we treat the inhomogeneous Neumann boundary value problem (|1.8I) . We 
proceed analogously to the Dirichlet case, that is, we hrst consider the linearized equation 
and then construct a solution of the nonlinear problem (11.81) by means of the implicit function 
theorem. 

Note that, in the Dirichlet case, the fact that the operator —A^: T>{A^) defined 

by (|3.3p has a strictly negative spectral bound iLemma 13.21) was crucial in order to show that 
the linearized equation (13.ip admits maximal regularity on M_|_, see the proof of Theorem 
13.31 In the Neumann case, due to the zero eigenvalue of —Aw: 'D{An) Lp{fl) with 

V{A]\f) = {u G Wp{Q) : dyU = 0 on T}, we cannot expect to obtain maximal regularity on M+. 
For this reason we consider —A7v,o: ^(Aw,o) ^ Lpfl{Ul), where D(A7v,o) = 7^(Aw) nLp^o(^^)- 
The spectrum of — A7v,o is contained in (0, 00 ), therefore we can prove maximal regularity of 
the homogeneous linear Neumann boundary problem on M_|_ analogously to the Dirichlet case. 
However, if we restrict ourselves to finite time intervals J = (0,T), then we do not necessarily 
need to use the realization —A7v,o in Lpfl{Q). In case of finite time intervals we use —Aw. As 
a consequence, we will prove global well-posedness of ()1.8p only if the data uq, ui, U 2 and g, 
h have zero mean whereas local well-posedness holds also for data with non-zero mean. 

4.1. Maximal Lp-regularity for the linearized equation. As in Section 13.11 let J = 

(0,T) or J = M+ and assume p G (1, 00 ). Here, for we f € Lp{J x D) we consider 

(oA — dt){utt — bAut — (?Au) = f in J x D, 

< {dyU, dyAu) = {g, h) on J x T, 

{u,ut,uu) = iuo,ui,U 2 ) on {f = 0} X D, 


(4.1) 
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where uo,ui,U 2 - —)■ M and g,h: J xT ^ M. are the given initial and boundary data, 
respectively. Analogously to the Dirichlet case we first represent (jd.ljl with g = h = 0 as an 
abstract evolution equation of the form 

(4.2) 

by setting 


dtv^+ =F, v^\0)=v, 


\N„.N 


N / 


..N 


V^ = 


u 

Ut 


- 

Un — 


Uq 

Ul 


and F = 


yU2 — c^A]\[Uo — bA]\[Ui 


yUtt - c^Anu - bAisfUt^ 
introducing the Banach space 

= P((Ajv)2) X V{An) X Lp{n) 
and defining the coefficient operator : F{A^) —)• via 



/ 0 -/ 0 \ 

(4.3) A^ = -c^An -bAN -I , V{A^) = V{{An)^) x V{{Am)^) x V{An). 

y 0 0 —aAw/ 


On one hand, in the following we will show maximal regularity of A^ on finite time intervals. 
On the other hand, as already pointed out, we are going to use the realization — Ajv,o of the 
homogeneous Neumann Laplacian. For this reason we introduce the Banach space 

X^’^ = V((Ajv,of) X V(Ajv,o) X Lp,o(^^), 

and the densely defined operator A^’^ : ^(A^’^) —?> X^’^, where Ajy has to be replaced by 
Aat^o in (SSI). 


Proposition 4.1. Let p € (1, cx)). There exists some v > 0 such that the operators v + A^ 
and V + A^'^ admit maximal regularity on M+. 


Proof. The result can be proved similarly to Proposition 13.11 For some a > 0 consider 


/al 

-I 




/ -al 

0 


0 

al — bAj\f 

-I 

and 

II 

-c^An 

—al 

0 

VO 

0 

al — a An I 



V 0 

0 

-all 


Clearly, the operator : X^ —>• X^ is bounded. Moreover, A^ : 'D(A)^) —)• X^ has max¬ 
imal Lp-regularity on M_|_ which is seen as in the proof of Proposition 13.11 by considering 
vt + A^v = F, no = 0 for n = (ni, ^2, t’s)"'" and F = (/i, /2, /s)'''. Explicitly, we have 

dtvi + avi - V2 = fi, ni(0)=0, 

dtV2 + aV2 - bANV2 - V3 = f2, V2{0)=0, 

dtV3 + av3 - a Awns = /s, ^ 3 ( 0 ) = 0 . 

Let F G Lp(M+; A-^). Now one solves stepwise the equations above, starting with the last 
one, to get a unique solution V 3 G VFp (1R+; Lp(fl)) n Lp(M+; T>(AAr)). Then f 2 + V 3 G 
Lp(]R+; T>(AAr)). Here, we need to employ Lemma [C.5I in order to obtain a unique solu¬ 
tion V 2 G VFp (M+; F(Ajv)) n Lp(M+; ^((Atv)^)). As in the Dirichlet case, the first equa¬ 
tion gives us a unique solution ni G Wp (M+; ^((Atv)^)). Altogether, since the condition 
F G Lp(]R+; A-^) implies existence of a unique solution v G lFp(M+; A'^) D Lp(M+; ©(A-^)) 
we conclude that A^ : 'D(A'^) —>■ X^ admits maximal Lp-regularity on M_|_. Finally, as in the 
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proof of Proposition 13.11 a perturbation argument implies that there exists some u > 0 such 
that V + + A 2 = v + A^ has the property of maximal Lp-regularity on . 

Maximal Lp-regularity of u + follows analogously by considering the operators ^4^’° 
and A^’^ which are equal to A^ and A 2 upon replacement of Atv by Aat^q and proceeding 
as above. □ 


Since u + ^4^^ has maximal Lp-regularity on M_|_, multiplication of (14.2j) with e shows 
that ^4^^ has maximal Lp-regularity on bounded intervals J = (0,r). 

Theorem 4.2. Suppose J = (0, T) is a finite time interval and let p € (l,oo). Then A^ : 
T>{A^) —>■ has maximal Lp-regularity on J and therefore 

{dt + A^,^,):W^iJ;X^)nLpiJ-ViA^))^LpiJ-,X^) x (X^,ViA^))i-i/p,p 

is an isomorphism. 

Next, observe that —A^’^ has a strictly negative spectral bound. This can be shown 
likewise to Lemma 13.21 since the spectrum of the negative Neumann Laplacian in Lp^o(^) is 
contained in (0, 00) and consists only of eigenvalues of hnite multiplicity. 

Lemma 4.3. The spectral bound of —A^’^ is given by s(— ^4'^’'^) = —oJq, where oJq = 
minjaA^, 6A^/2, c^/6}. Here, denotes the smallest non-zero eigenvalue of —An^q. In 
particular, i/Re(A) < ojq , then A G p{A^’^). 

By means of Lemma [4. 3 1 one shows that A^’^ has maximal Lp-regularity on M+. For details 
we refer to the proof of Theorem 13.31 

Theorem 4.4. Let p G (l,oo) and uj G Then A^’^ : 'D{A^’^) —)> X^’^ has maximal 

Lp-regularity on M_|_ and therefore 

{dt + A^’^,-ft): e-^{W^{J;X^’^)nLp{J-,V{A^’^))) 

^e-Lp(J;X^’0) X {X^^^,V{A^’^)fi_yp^p 

is an isomorphism. 

Theorems 14.21 and 14.41 immediately yield optimal regularity for (j4.ip with homogeneous 
boundary conditions, i. e. <7 = h = 0. 


Corollary 4.5. Let p G (l,oo)\{3} and consider the homogeneous Neumann boundary value 
problem 


(4.4) 


(aA - dt){utt - bAut - c^Au) 
< {dyU, dyAu) 

{u,ut,utt) 


f in J X Ll, 

(0,0) on J X T, 

{uo,ui,U2) on {t = 0} X Ll. 


(i) If J = (0, T) is finite, then (14.4p admits optimal regularity in the sense that there exists 
a unique solution 


u G E,(J), E„(J) = W^{J- Lp(54)) n W^{J- W^{n)) 

if and only if f ^ Lp{J x Al), uq G Wp{Q), ui G Wp U 2 G Wp and the 

initial and boundary data are compatible, that is, we have dyUofi = dyAuo\r = dyufir = 
0 and, if p > 3, also dyAui\Y = dyU 2 \T = 0 in the sense of traces. 
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(ii) If J = M_|_, then for every ui € (0 ,lv^) with co^ = mm{aX^,bXf /2,cP/b} we have that 
dMl) admits optimal regularity in the sense that there exists a unique solution 

u G e-‘"E„,o, E .,0 = Wp3(M+; Lp,o(f^)) n Ty;(M+; W^{n) n Lp,o(^^)) 

if and only if f G e“'^Lp(M+; Lp,o(^^)), uo G Wp{Q)r\Lpfi{Q,), m G Wp~^''^{Il)r\Lpfi{Il), 
U 2 G Wp n Lpfi{Il) and the initial and boundary data are compatible. 

Proof. Assertion (i) follows immediately from Theorem 14.21 Analogously to the proof of 
Corollary 13.41 one verifies that G Wp{J; X^) n Lp{J]T>{A^)), F G Lp{J;X^) and Vq G 
[X^,'D{A^))i_i/pp imply u G Eu(J), / G Lp{J x 14) and the desired regularity of the initial 
values, respectively. Based on Theorem 14.41 the second claim follows analogously. □ 


Finally we arrive at our global optimal regularity result for gT]). As in the Dirichlet 
case, sufficiency is shown by a combination of an optimal regularity result for the linearized 
Westervelt equation and a higher regularity result for the heat equation. 


Proposition 4.6. Let p G (l,oo) \ {3} and define ujq = min{aA]^, 6A]^/2, c^/6}. Then for 
every to G (0,a;o) the linear initial boundary value problem 


(4.5) 


' (aA - dt){uu - bAut - c^Au) 
< {dyU, dyAu) 

{u,ut,utt) 


f in M_|_ X n, 

{g, h) on R+ x T, 

(uo,ui,U2) on {t = 0} X Ll. 


has a unique solution of the form u{t,x) = v{t,x) + w{t), where 
V G e-‘"E„,o, E,,o = bFp3(M+; Lp,o(f^)) C Ty;(M+; W^{n) n Lp,o(bl)), dfw G e-‘"Lp(M+) 


if and only if the data satisfy the conditions 
(i) / G e-“Lp(M+ X n), 

(ii) uo G W^{n), ui G bFp"^/^(14), U 2 G Wp~^^^(n), 

(iii) g G e-Fp,,, Lp{T)) n bF;(M+, Wp""^/'’(r)), 

h G e—F;,,, = fFp^/'-^/'^(]R+; Lp(r)) n fFpHM+; fFp^“^/'’(r)), 

(iv) dyUolr = g\t=o, dyAuo\r = h\t=o, dyUi\T = gt\t=o and, if p > 3, also 5j.Aui|r = ht\t=o, 
dyU 2 \T = gtt\t=o in the sense of traces. 

Moreover, the solution fulfills the estimate 

ll'*^lle-“E„ ~ WfWe-'^Lp + Ilfl'lle-^F,,,^ + ll^lle-‘^Fh,„ + ll'*^o||iy4 + ||ui||^4-2/p + ||t(21|^^2-2/p. 


Proof. It is not surprising that the proof of necessity can be done similarly to Proposition 
ESI Assume that u(t, x) = v{t, x) + w(t, x) is a solution of (14.5p with u G e ‘^Eu and 
dfw G e“^Lp(M+). Since, apart from having zero mean, v has the same regularity as u in 
Proposition 13.51 we are conclude that e“*/ = —e^*(u+ u;)t4t — a6A^(e‘^*ut) — ac^A(e‘^^u) + (a + 
b)A{e‘^*vu) + c^A^i&’^vt) G Tp(M+ x 14) and (i) is readily checked. Moreover, w, wt and wu 
are just time-dependent and thus constant at t = 0, hence the regularity of the initial values 
(ii) can be shown as in the Dirichlet case. Moreover, the regularity of the boundary data (iii) 
is obtained from the spatial trace (12.2p with the same choices of k and I as in the proof of 
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Proposition 13.51 and setting j_b = 1. Concerning (iv) it is straightforward to show 
d,yUo\r = g\t=o in 5i.Auo|r = h\t=o in 

= gt\t=o in d^/\ui\Y = ht\t=Q in WI~^/p{T) if p > 3, 

dvU2\r = gtt\t=o in Wp~^^'P{T) if p > 3. 

Next, we show that conditions (i)-(iv) are sufficient for the existence of a unique solution 
u{t,x) = v{t,x) + w{t) of (j4.5j) such that v € and wut G e“‘^Lp(M+). As in the 

Dirichlet case, we interchange the order of differentiation on the left-hand side and consider 
the subproblems 



iftt - bAipt -c^Aip = f 


in M_|_ X H, 

(4.6) 

duP = ah- gt 


on M_|_ X F, 


= (aAuo- 

ui,aAui — U 2 ) 

on {f = 0} X H, 

and 

f aAu — ut = (p 

in M_|_ X H, 


(4.7) 

J d^u = g 

on M_|_ X F, 



u = uq on {t = 0} X fl. 


From condition (ii) we obtain uAuq — ui G PPp(fl) and aAui — U 2 G Wp Furthermore, 

(iii) implies ah — gt G On the strength of Lemma 12.71 we obtain that (|4.6p admits 

a unique solution of the form ^p{t,x) = (pi{t,x) -|- (p 2 {t) with ipi G e“^Wu,o and £ 

e“‘^Lp(M+). We now make the ansatz u{x,t) = v{x,t) -|- w{t) such that v{-,t) = 0. Applying 
/o to aAu — ut = we deduce that w solves the ordinary differential equation wt = 
—(/?2 + with w{0) = Uq. Hence dfw G e“‘^Lp(M+). Moreover, u is a solution of 

(4.8) aAu — vt = ipi + a|r||fl|“^^ in 14, d,yV = g onT, u(0) = uq — uq. 

In order to apply Corollary IC.41 we first note that the right-hand side (/?i + a|r||14|“^^ belongs 
to e~‘^Wu since g only depends on time and belongs to (M+). The rescaled function 

Va{t,x) = av{t/a,x) should solve the system 

(4.9) Ava - dtVa = ^^1 + a|r||H|“^5 in H, dyVa = ag on F, Va(0) = auo - auQ. 

Hence the compatibility condition ()C.7p becomes 


— / ((^1-|-o|F||H| g) dx + / agdS = 0 


and is clearly satished. Therefore Corollary 1C.41 yields a unique solution Va G e“^E„^o of 
problem (|4.9I) and thus u = v + w solves problem (|4.5I) . 

Finally, uniqueness follows by considering two solutions of (14.51) . the difference u of which 
solves (j4.4l) with / = 0 and uq = ui = U 2 = 0, hence a = 0 and the proof is complete. □ 


Proposition 4.7. Let p G (l,oo)\{3} and let J = (0,T) be a finite interval. Then the linear 
initial boundary value problem has a unique solution 

u G Wp3( J; Lp(L!)) n Wp^J; 

if and only if conditions (i)-(iv) from Provosition 14.61 hold with replaced by J. 
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Proof. Uniqueness and necessity are shown likewise to Proposition 14.61 For the proof of 
sufficiency one considers (j4.6ji and ()4.7ji on (0,T) instead of M+. Remark 12.81 and Corollary 
1C.41 then imply existence of a unique solution u G Wp{J] Lp{n)) n VFp (J; Wp{PL)). □ 

4.2. Local well-posedness, global well-posedness and exponential stability. We now 

arrive at our well-posedness results for the Neumann problem (II.8p . If we allow the (suffi¬ 
ciently small) initial and boundary data to have non-zero mean, we are only able to prove 
well-posedness of (jl.8p on finite time intervals J = (0,T). On the other hand, if the (suffi¬ 
ciently small) data have zero mean, that is, ho = hi = 1x2 = 0 and 5 = /i = 0, then we obtain 
a globally well-posed solution which decays exponentially fast to zero. 

Theorem 4.8 (Local well-posedness). Let J = (0,T) for some T < 00 and p > max{n/4 -|- 
1/2, n/3}, p 7^ 3. Suppose 

uo € W^{n), xxi G U 2 G W2-2/p(L!), 

(4.10) g G e-‘"F,,,(J), F,,,(J) = {J, Lp{T)) C 

h G e-‘"F;,,,(J), F,,,,(J) = Lp{T)) n W^iJ- 

such that duUo\r = 9\t=o, ^lyAxxoIr = h\t=o, dyUi\T = gt\t=o and, if p > 3, also duAuilr = 
ht\t=o, d,yU 2 \r = gtt\t=o in the sense of traces. 

There exists some p > 0 such that if 

IIS'llFg,,. + + IkollvE/ + lkl|liy4-2/p -I- ||xX2||^2-2/p < p, 

the nonlinear Neumann boundary value problem (USD admits a unique solution 

(4.11) XX G E^J) = iy3(J; Lp{n)) O W^{J; w^{n)) 

which depends analytically (in particular continuously) on the data (|4.10l) with respect to the 
corresponding topologies. Moreover, conditions (|4.10l) are necessary for the regularity of the 
solution given in (14.lip . 

Theorem 4.9 (Global well-posedness - the Neumann case). Let p > max{n/4 -|- 1/2, n/3}, 
p 7^ 3 and define loq = min{aA)^, bX^/2, (?/b). Suppose 

no G w^{a) n Lp,o(f^), xxi G n Lp,o(bi), XX 2 G n Lp,o(Ll), 

(4.12) g G e-‘"Fp,,,o, Eg,,,o = Lp,o(r)) C W^iR+, C Lp,o(r)), 

h G e--F,,,,,o, Lp,o(r)) n WpHM+; W^-^ppr) n Lp,o(r)), 

such that duUo\r = g\t=o, ^i^AxxoIr = h|t=o, dj,xxi|r = 5t|x=o and, if p > 3, also d^^Axxilr = 
^t|i=0; 5i/XX2|r = gtt\t=o in the sense of traces. 

Then for every oj G (0,a;o) there exists some p > 0 such that if 

ll5'lle-“Fg,„,o + ll^lle-“Fh,„,o + ll^ollvE/ + ll^4-2/p -|- ||n2||^2-2/p < p, 

then the nonlinear Neumann boundary value problem (USD admits a unique solution 

(4.13) n G e--E,,o, E„,o = (E+; Lp,o(Ll)) 0 lUpHM+; W^{a) n Lp,o(0)), 

which depends analytically on the data p4.12p with respect to the corresponding topologies. 
Moreover, conditions (14.121) are necessary for the regularity of the solution given in (14.131) . 
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Remark 4.10. In case of Theorem 14.81 we define rt* to be the solution according to Propo¬ 
sition [Q which satisfies ()4.ip for the data (/ = 0, g, h,uo,ui,U 2 ) und suppose u, satisfies 
homogeneous boundary and initial conditions. The solution is then of the form rt = w* + u, 
and u, is found by the implicit function theorem. The claim then follows likewise to the proof 
of Theorem 13.61 

However, if we want to prove global well-posedness, we need to use Proposition 14.61 for 
the linearized equation, where for given data (/ = 0, g, h,uo,ui,U 2 ) according to (ii)-(iv) the 
solution is of the form u{t, x) = v{t, x) + w{t), where v G has zero mean and w is only 

time-dependent. If re 7 ^ 0, the term in the nonlinear right-hand side of (II. 8 p causes 

problems. Recall (13.121) and note that due to Proposition |4T] we in fact have tt* = with 

u* G e“‘^Eu^o and G e“‘^Lp(M+). Then ui*, dtWi, and are in general not contained 
in e“‘^Lp(M+) and thus (I3.13p fails. However, if we assume that the data {g, h, uq, tti, U 2 ) have 
zero mean then rc* = 0 and, since tt* = u* in this case. Theorem follows analogously to 
the result on global well-posedness for the Dirichlet boundary value problem (Theorem 13.61) . 

Finally, provided the data have zero mean, we obtain the following result on exponential 
stability for the Neumann problem (11.81) . 

Theorem 4.11 (Exponential stability - the Neumann case). Under the same assumptions as 
in Theorem KM the solution u decays exponentially fast to zero as t —>■ 00 in the sense that 

\Wi't)\\wf + \\ut{t)\\^^^- 2 /p + \\uu{t)\\■^y 2 - 2 /p < Ce t > 0, 

for some (7 > 0 depending on the boundary and initial data g, h, uq, ui and U2- 

Proof. Note that we have u G e^^E^^o ^ e“‘^E„, therefore the result follows likewise to 
Theorem 13.71 □ 

Appendix A. The Neumann Laplace operator 

Let p G ( 1 , 00 ) and assume, as always, that Ut C M"" is a bounded domain with smooth 
boundary P = dit. The homogeneous Neumann-Laplacian is given by 

—Atv: Lp{Ut), 

u i-A —Au, 

where 'D{Aj\f) = {u € Wp{Ut)\ dyU = 0 on T}. It is well-known that — Atv has compact 
resolvent and that its spectrum ct(—Atv) is a discrete subset of [ 0 , 00 ) consisting only of 
eigenvalues {Xn)n>o with finite multiplicity. In particular, 0 = Xq G it(— A^r) is an isolated 
eigenvalue of — A^v. We seek for a realization of the Laplace operator with homogeneous 
Neumann boundary conditions such that the spectrum is contained in [X^, 00 ) where X^ > 0 
is the smallest non-zero eigenvalue of — Atv. 

In order to remove the zero eigenvalue we will use several results from Appendix A in 
|Lun95] . In what follows. A: T>{A) C A —> A denotes a linear closed linear operator whose 
domain T{A) is dense in the real or complex Banach space A 7 ^ {0}. We say that a subset 
(Ti C (y{A) is a spectral set if both, ai and (t{A) \ ui are closed in C. Let ui be a bounded 
spectral set and let 02 = <t(A) \ ui. Since dist((7i, < 72 ) > 0, there exists a bounded open set 
O such that ai <Z O and O n cr 2 = 0. We may assume that the boundary 7 of O consists 
of a finite number of rectifiable closed Jordan curves, oriented counterclockwise and define a 
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linear bounded operator P by 

The following result shows how find a realization A2 of A such that cr{A2) = cr{A) \ ci. 

Proposition A.l ( |Lun95l Proposition A.1.2]). Let ai be a bounded spectral set. Then the 
operator P is a projection and P{X) is contained in 'D{A'^) for every n G N. Moreover, if we 
set Xi = P{X), X2 = (/ — P){X) and define the operators 

Ai: Xi —>■ Xi, tt i-A Au and A2 : 'D{A2) = P{A) D X2 ^ X2, n i-)- Au 

then 

(t{Ai) = (Ti and cr{A2) = ( 72 - 

The crucial point is thus to determine the space X 2 . In case ai = {Aq} where Aq is an 
isolated point of (t(A) and a pole of R{-, A) the following result helps to determine the spaces 
Xi and X 2 . 

Proposition A .2 f |Lun95[ Proposition A.2.2 and Corollary A.2.4]). If Xq is an isolated point 
of a (A) and a pole ofR{-,A), then the following are equivalent: 

(i) Xi = N{XoI - A) 

(ii) X 2 = i?(Ao/ - A) 

(hi) Aq is a simple pole of X R{X, A) 

(iv) R{XqI — A) is closed and X = N{XoI — A) © R{Xo — A) 

(v) N{XoI - A) = N{XoI - A)2 

We now apply the foregoing results to the strong Neumann-Laplacian —Atv and set Aq = 
Xq = 0, hence ai = {0} and 1 T 2 = a{—Aj\f) \ {0}. Then we clearly have a 2 C [A]^, 00 ), where 
X^ is the smallest non-zero eigenvalue of —Aat. 

Lemma A. 3 . The spectrum of —A^ consists only of poles of X i?(A,—A^v). 

Proof. Since —A^r is closed, densely defined and has compact resolvent, the result is an 
immediate consequence of Corollary IV. 1.19 in |ENf)0] . □ 

We introduce the space = {u G Lp{Il): u is constant} and start with the following 
observation. 

Lemma A. 4 . We have N{A]\[) = 

Proof. Let u G X(A 7 v), i-e. u G Wp{Ll), dyU = 0 on T and —An = 0. For sufficiently large 
fj, > 0 the map p — An : I1(A})^^) —>■ I1(A-^) is an ismorphism vor every j G Nq. We write 
0 = Anu = piu — {pL — An)u and obtain u = {p — An)~^pu. Hence, if n G V{An) then 
we have u G 'D{{Aj\f)‘^) and altogether conclude u G V{{An)°°). The Sobolev embedding 
Wp{Ll) ^ IT|(H) holds for sufficiently large k. Therefore we in fact have u G IT|(H) and 
calculate 

0 = — / Auudx= / Xu-Vudx— / dyuudS = \\V, 

Jo. Jq Jt 

hence u G Conversely, every function u G trivially satisfies —An = 0. □ 

Lemma A. 5 . We have X(Ajv) = N{{An)‘^)■ 
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Proof. Let u G i.e. u € and {—Af^u = 0. Then 

0 = / {—A)‘^uudx = — / V{Au)-Vudx+ d^{Au)udS 
Jn Jn Jr 

= / AuAudx = WAuWl^, 

Jn 

hence —Au = 0 and thus u G N{Aj^). Conversely, let u G N{Aj\f). Then (—A)^u = 0 and 
—Art G I?(Ajv)- n 

Before we proceed, let us recall the space Lp^o(^) = {rt G Lp{Q): JqU = 0}. 

Lemma A.6. ITe have Lp(fl) = © Lp^o(Q) as a topological direct sum. 

Proof. We consider the map P: Lp{Q) —)• Lpfi{^), u u — {u)q where {u)q = J^udx. 
It is straightforward to verify that for rt G Lp{Q) we indeed have j^u— {u)^dx = Q and hence 
Pu G Lpfl{Q). Furthermore, P^u = Pu— {Pu)q = Pu implies that P is a projection. Finally 
N{P) = since u — (u)^ = 0 if and only if u G Lp{Q) is constant. □ 

Therewith, by means of Proposition IA.2I have determined the space X 2 = which 

gives us the following positive realization of the Neumann Laplacian. 

Theorem A.7. The spectrum of the closed and densely defined operator 

u !->• —Au, 

with P(Ajv,o) = n Lpfi{Q) is a discrete subset of [A^,oo), where > 0 is the 

smallest non-zero eigenvalue of—A^. Moreover, (t(—Aat^q) consists only of eigenvalues with 
finite algebraic multiplicity. 

Appendix B. Traces and mixed derivatives 

In this section we consider the temporal trace operator in some anisotropic fractional 
Sobolev spaces. Furthermore, we present some mixed derivative embeddings for such spaces 
which are needed for proving suitable mapping properties of differential operators. 

A bounded linear operator r: X ^ Y between Banach spaces X and Y is called a retraction, 
if there is a bounded linear map r^:Y —>■ A such that rr^ = ly. Thus r is surjective and r^ 
is a bounded right-inverse for r. The map is called a co-retraction for r. 

The following trace theorem can be derived from [DB84t Lemma 11], [Lunhb] Section 2.2.1], 
|Ama951 Proposition III.4.10.3]. 

Theorem B.l. Let A be the generator of a bounded analytic semigroup in- 0 . Banach 

space X such that A: T>{A) —>• X has a bounded inverse, let p G (l,oo) and let T>A{oi,p) ■= 
{X,V{A))a,p for a G (0,1) and Va{1,p) := V{A). 

Then, for every a G (1/p, 1], the trace operator 

7t = -14=0: n Lp(M+;PA(a,p)) - l/p,p) 

is a retraction, the operator 

Ra: uo ^ {t e"*^uo) , - l/p,p) ^ 1T"(M+; A) n Lp(M+;P^(a,p)) 

is a co-retracton for cmd the following embedding is continuous. 

1T“(M+; A) nLp(M+;pA(a,p)) -A BUC{R+-,VAia - l/p,p)). 
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This theorem can be applied to the spaces 1T"(]R+; Wp{Q] E)) n Lp(]R_|_; E) for 

a G ( 1 /p, 1 ], s G [ 0 , oo), provided that in the case a < 1 the number s + 2 q: is not an integer 
and provided that s + 2a < 2k. Here and in the following, we assume that E is a Banach space 
of class 'H'T and has property (a), where we refer to |KW04| and [KS12] for the definitions of 
such spaces and additional information. For instance, any Hilbert space is of class TiT with 
property (a) and these properties are inherited to closed subspaces and isomorphic spaces. 
Moreover, the space Lq{Q,A,fJ.]E) on a a-finite measure space with q G (l,oo) is 

of class T-LT and has property (a). For s G (0, oo) and p, q € (l,oo), the Sobolev-Slobodeckh 
spaces VFp (M"'; Fi)), the Bessel potential spaces E[p{W^-,E) and the Besov spaces 
and are also of class T-LT with property (a). 

Let us indicate how Theorem IB . 1 1 can be applied. Let first be an extension operator from 
H to M"' which acts as a bounded linear operator Wp{Q]E) Wp(ME-,E) for all t G [0, 2A:]. 
Such extension operators are defined in |AF03| for t G No and their boundedness for t ^ Nq 
follows from real interpolation. Then it remains to study Squ in M_|_ x M”. In this situation 
the operator A = l — A in X = ITp(R"';£') with domain ITp+^(]R"'; Fi) has the required 
properties. Indeed, [DHPD^t Section 5] covers the case s = 0 and an abstract result of Dore 
[Dor99| covers the case s G (0,2) \ {1}. The remaining cases follow by means of isomorphic 
mappings, interpolation and taking fractional powers. Hence 'DA{a,p) = VFp'’'^“(M"'; Fi) and 

T>A{ci — l/p,p) = ^^^(M"';Fi). Then the temporal trace operator can be rewritten as 

7 i: n i-A ((Ton)|t=o)|o and acts as a bounded linear operator 

IT;(M+; w;{n- E)) n Fp(M+; ^ ^.+ 2 a- 2 /p(^. 

For the boundary spaces HT*(M+; hFp (F; Fi)) n Lp(M+; ITp+^"(r; Fi)) we use a common 
retraction r: Wp{M.'^~^]E)^ — Wp{T;E) for all t G [0,2F] with some A" G N. A co-retraction 
for r can be constructed by means of a partition of unity for F and local parametrizations of 
F over subsets of as in the proof of Lemma IC.61 Then the temporal trace operator can 
be rewritten as 74 : u i->- r{{r^u)\t=o) and maps 

W^iR+-,W;{T-,E))nLp{R+-W;+^^{T-,E)) ^ E). 

In order to construct functions with prescribed initial values, we consider an operator 
A: T>{A) C A —>■ A as in Theorem IB . 1 1 and define the spaces 

VA{k + a,p) := A~^VA{a,p) = {V{A^),V{A^"^^))a,p for F G No, a G [0,1], p G (l,oo). 
Then Theorem IB . 1 1 and the identity dye~^^ = —Ae~^^ = e~^^A yield the following result. 
Corollary B. 2 . Let k G Nq, a G (1/p, 1] and p G (1, 00). Then the operator 

Ra : u i-a (t i-a e~^^u) , VA{k + a — l/p,p) ^ ITp’''"(M+; A) n Fp(M+; Da(F + a,p)) 
is a bounded right-inverse for 7 ^. 

We next deal with higher order initial conditions. 

Lemma B. 3 . Let 7 / := (9j-)|t=o and let I G Nq, m G N with m > I 1. Then the operator 

(70,7/,..., 7^): H 7 (M+; A) n Lp{R+-V{A^)) ^ J] ■^o^A(m - j - l/p,p) 


is a retraction. 
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Proof. For j G {0,1,..., 1} and x € X we define 

:= A~^X for f > 0 , 

■’ ^ —^ 1=0 

where, for each j, the I + 1 numbers Cij (z G { 0 ,1 ,... , /}) solve the linear system 

oCd(~(^ + *))”" = ^rnj for m G {0,1,... , /}. 

By using Vandermonde’s matrix 

/I 1 ... 1 \ 


V = 


1 + k 


with 


fe(fc+l) *1 —r 

detV = (-l)—nj!/0, 
i=i 


2^ ... {l + kf) 
the numbers Cij are given by (coj ,... ,cij)^ '■= Hence 

{df^Sfx){0) = = 6mjX for m G {0,1,... ,1}. 


From Corollary IB.21 we infer that acts as a bounded linear operator 


Sf : VA{m - j - l/p,p) ^ H7(M+; V) n Lp 
Therefore the desired co-retraction is given by 

S^{xo,xi,... ,xi) := V]. S\ 

^ — ^j=0 ■ 


; 'D(A^)) for m G N, m > J + 1. 


Xj. 

=n J 


□ 


Theorem B .4 (Mixed derivative embeddings). Let n G N, p G (l,oo), G [0,oo), T,a ^ 
(0, oo), 6 € (0,1), S G LIT, J = M. or J = {0,T) for T G (0, oo] and let H be the whole space 
M” or a bounded domain with smooth boundary or a compact smooth hypersurface o/M”. 
Then the following embeddings are continuous. 

(B.la) J; H;{n-E)) n i/*(J; E)) ^ {n-E)), 

(B.ib) BY;{J;H;{n;E))nH^p{j;B;y{n-,E)) ^ BY;^j-,H;+^^-^>{n-,E)), 

(B.lc) BY;iJ- E)) n J; E)) ^ J; E)). 

Proof. We adapt the proof of |MS12l Proposition 3.2], It is sufficient to consider the case 
J X n = R X R” since the other spaces are retracts of corresponding spaces over R x R"'. 

In the ground space X = Hp{Hf) := Hp{R; E)) we consider the operators 

H = (I - V{A) = F*+"(R; E)), 

B = {1- V{B) = Hl,{^-H;+%mA-E)). 

Here A„ denotes the Laplacian in R"" with respect to the spatial variable. It follows from 
|DHP03[ Theorem 5.5] (see also |MS12[ Lemma 3.1]) that the operator 

Jr := (I - A ™)^/2 . h;+p{W^-E) ^ f;(r™; F) 

is invertible and has a bounded functional calculus and thus bounded imaginary pow¬ 
ers in Hp(R™;T) for all m G {l,n}, r,p € [0, oo) and all Banach spaces F of class LIT. 
The latter property implies that its fractional powers {9 G (0,1)) have the domains 
V{jf) = [H'^{W^-F),hYP{W^] F)]e = F) by |DHP03[ Theorem 2.5] and complex 

interpolation. 
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By choosing {m,r,p,F) = E)) or {m,r,p,F) = {n, s,a, E)) and em¬ 

ploying Fubini’s theorem, we see that A: F^A) —)• X and B: F{B) —)• X are invertible and 
have bounded 71°° functional calculi. Sobolevskii’s mixed derivative theorem |Sob75l Theorem 
6] implies that F{A) D F{B) 'D{A^B^~^) n V{B^~^A^) which proves (IB.lap . 

For proving (jB.lbp we apply (IB.lap and real interpolation to the space 

for sufficiently small e > 0. For t± := (1 ± e)T and a± := (1 ± e)a and 9± G (0,1), we obtain 


We choose 9± such that (1 — 9±)a± = (1 — 9)a, that is, 9± := {9 ± e)/(l ± e). Then 
9±e± = 9t ±€T and it remains to apply the real interpolation functor (•, •)i/ 2 ,p to 

Indeed, interpolation of the right-hand side yields 


(Z_, Z+)i/ 2 ^p ^ 


/ ut-\-6T—eT 

\ P 




)) 


1/2,P 




For an interpolation of the left-hand side we write = F{L) and Z_ = <^)/(i+'^)) 

where the operator L = (1 — clj 2 ^(i+e)T /2 -(- (i — is considered in the ground space 

Zq := Hp(Hp). Then the reiteration theorem (^ [Lun95l Remark 1.2.16]) yields 


(Z_,Z+)i/2,p = (P(l(i-^)/(i+^)),P(L))i/2,p 

= Vl (1/2 + (1 - e)/(2 + 26),p) = Bj,-]p^{H;) n 77*(B;+"). 

Hence (IB.lbp is proved. The proof of (IB.lcp is similar and therefore omitted. □ 


Appendix C. Higher regularity for the heat equation 


We study the regularity of solutions of the heat problem 

' {dt + Pb — A)tt = / in J X H, 

(C.l) JBU = g on J X F, 

u|i=o = uo in Q. 


Here J is a bounded interval (0, T) or the half line (0, oo) and 11 is a bounded domain in M"', 
n G N, with smooth boundary F. For B G {D,N}, let pB be a real number and let 

ID ■= -Ir, 7N ■■= (9i.-)lr = i' ■ (V-)|r, 7t := (^/•)k=o, 7t := 7? 


denote the Dirichlet, the Neumann, and the temporal trace operators, respectively. Again we 
let Xq > 0 denote the smallest eigenvalue of — Apj and > 0 denote the smallest non-zero 
eigenvalue of — Atv- We will prove the following regularity result. 


Theorem C.l. Let B G {D,N}, Jb = 0, Jn = 1, pd ^ (“-^^jOo), Pn € (0, oo), I G No, 
/c G N and p G (1, oo) such that Jb/'T- + 3/2p ^ 1. Then problem (jC.ip has a unique solution 


(C.2) 


u G := J; Tp(H)) n W^{J; W^\n)), 
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if and only if the data (/, 5,tto) satisfy the regularity conditions 


(C.3a) 

(C.3b) 

(C.3c) 


/ e n W'(J; -2(0)), 

g G j. j. ^^2fc-iB-l/p(r)), 

■uo e 74E = 


/■2k 

P 

/2k- 


Ty2fc-2/p(57) if 1 = 0^ 


and the compatibility conditions 


(C.4a) 

(C.4b) 


V + (A - //B)^^j-i e 


w^^{n) 

bbp 2 ('+^-^')- 2 /P( 0 ) 


Jt9 = iBUj for j € No, j <l + k- jB/2 - 3/2p. 


forj G Nn [1,/ - 1], 
for j G N n [Z, / + /c — 1], 


Remark C.2. (i) The space E*^7 ^ — 2 ) ig the standard parabolic solution space, 

(ii) If Z = 0 and J x O is the half space M+ x M"’ or the wedge M+ x M”, then E*^’^ is the 

anisotropic space {J x 0) with weight i' = (2,1,... , 1) in the sense of [Ama09] . 

This fact will be used in the construction of functions with prescribed boundary values, 
(hi) We exclude the case js/2 + 3/2p = 1 in order to avoid the more complicated trace 

spaces 71)14^3^^2 (^) 77 vl 4 ^ 3 ^^( 0 ). 

(iv) The additional regularity conditions (IC.4ap follow from the non-triangular structure of 
the space E^’^ in the case I > 1 and are derived in Subsection 1C.11 For j > l + l, 
formula (IC.4ap does not contain additional regularity conditions and should be merely 
understood as the definition of the functions Uj, which appear in ()C.4bp . 

(v) Every solution satisfies the higher order boundary conditions 

(C.5) 7 BA-^’+^n = {dt + HB)gj - f ='■ Qj+i for j G No n [0, k-2], with go := g. 


With the temporal trace theorem and Ui = 7^11 we obtain 

(C. 6 ) 'YB^^Ui = YiPj for j € No n [ 0 , /c - 1 ], i G No n [ 0 , Z + /c - j - jB/2 - 3/2p]. 

These equations are no additional regularity or compatibility conditions but follow from 
(IC.4p . (jC.5p . by induction over j G Nq. Indeed, suppose that = 7t9j for all i 

and some j. Then (IC.hp . the induction hypothesis, and (IC.4I1 yield 

7t5i+i = 'll {(.dt + hB)9j - iB^^f) = ll^^gj + 9Bll9j - iBll^^f 
= 7BA%i+i + HBlB^^Ui - iBll^^f 

= 7 BAf {-/If + (A - //B)'Wi) + tiBlB^^Ui - 'YBll^^f = 7sAf+^ni. 

(vi) In the case I = 1, k = 2, B = D, the compatibility conditions read as 

ui := f\t=o + (A - giD)uo G VFf (12), 

U 2 := dtf\t=o + (A - fiD)ui G Wf-2/P(57), 
g\t=o = ^^o|r, 

dtg\t=o = f\t=o,r + ((A - /iD)tto)|r, 

d^g\t=o = dtf\t=o,r + ((A - //D)'Wi)|r if P > 3/2. 
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(vii) The corresponding result for B = N has the compatibility conditions 

ui := f\t=o + (A - fiN)uo € Vhf (O), 

U2 ■■= dtf\t=o + (A - fiN)ui £ Wf-2/P(0), 

g\t=o = duUQ\r, 

dt9\t=o = di,f\t=o,r + 5y((A — /iAr)uo)|r, 

dig\t=o = d^dtf\t=o,r + ^^.((A - g.N)ui)\r if p > 3. 


Aiming at stability for Neumann boundary conditions, we will also prove the following 
result, where we consider the subspace Lpfi{Q) := {/ G Lp{Q) : f{x)dx = 0}. 

Corollary C. 3. Let B = N, pn £ (—A^,oo), I £ Nq, k £ N, p £ (l,oo), p ^ 3. Then 
problem (jC.ip has a unique solution 

u£E^^^ :=E^'’^nLp{J;Lpp{n)), 

if and only if the data {f,g,uo) satisfy the regularity conditions (IC.3P and the compatibility 
conditions ()C.4p and 

(C.7) / uo{x)dx = 0, / f{t,x)dx+ / g{t,x) dS{x) = 0 for t £ J. 

Jn Jn Jr 


Next, we study the original heat problem 

^ {dt - A)u = f 

(C. 8 ) 


IBU = g 
u\t=o = ^^0 


in X Ll, 
on R+ X r, 
in Ll. 


Corollary C.4. (i) Let B £ {D,N}, po £ (—A^,oo), pjy £ (0, oo), I £ No, k £ N, 
p £ (l,oo) such that js/‘^ + 3/2p 1. Then problem (|C. 8 h has a unique solution u £ 

if and only if the data {f,g,uo) satisfy the regularity conditions 


(C.9) (/, g, uo) £ x e^^ysE^’^ x 

and the compatibility conditions 




ifl>f, 
if 1 = 0, 


(C.lOa) 

(C.lOb) 


'“j := It V + Auj-i £ 


r ITf(O) /ori€Nn[l,/-l], 

\ for j £Nr\[l,l + k-l], 


Jt9 = iBUj for j G No n [0, / + A; - Jb/^ - 3/2p]. 


(ii) Let B = N, p^ £ (—A^,oo). Then problem (IC.SP has a unique solution u £ c^^Eq^ 
if and only if the data (/, 5 ,uo) satisfy the regularity conditions (IC.9P and the compatibility 
conditions (jC.lOp . (IC.7I) . 


Proof. In problem (jC.ip we multiply /, g with e^®*, so that 

This shows that solves ()C. 8 p for (e^®*/, uq) if and only if u solves the shifted 

problem (jC.ip for {f,g,uo). Hence Theorem 1C.II and Corollary 1C.31 yield the assertions. □ 
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C.l. Compatibility conditions. By means of the results from Section [B] it is not difficult 
to verify that the regularity conditions ()C.3() are indeed necessary for u G E^’^. Let us now 
derive the remaining compatibility conditions. 

First, any function u G E^’^ satishes the initial regularity conditions 


(C.ll) 




foriGNon[0,Z-l], 

^ 2 l+ 2 k- 23 - 2 /p^^^ for j enon[l,l + k-i]. 

-1 


If u solves (jC.ip with data {f,g,uo), then an application of ^ to the heat equation yields 


dtu = ^t f + {^-k-B)di 


-1 


In particular, the initial values Uj := dlu\t=Q are given in terms of / and uq by 








Itf + (A - gsY Uq for j G N n [1, / + /c - 1]. 


Then (IC.llI) implies that the data (/, uq) and uj must satisfy 

rwf(L!) for jGNn[l,/-l], 

\ for j G N n [/, / + A: - 1]. 

For j G [1,Z] this is indeed an additional condition, since / merely satisfies 


'“j = 7t V + (A - G 




W, 


2k-2 


(L!) 


y^2l+2k-2j-2/p( 


llg = iBUj G 


for j G No n [0,1], 

\n) for j G Nn [/ + 1,Z + A:-1]. 

The conditions for 1 + 1 < j < I + k — 1 then follow from the regularity of f, uq, ..., Ui 
and could therefore be omitted in (|C.4ah . but we keep them there as a dehnition of u^+i, 
..., ui^k-i- Indeed, these functions still admit traces on F. By differentiating the boundary 
condition ygu = g with respect to time, we obtain 

r Wf-^'^-i/PpF) for J G No n [0, / - 1], 

\ w2(^+^-j)-Js-3/P(r) for J G No n [Z, / + A: - jB/2 - 3/2p]. 

This shows that ()C.4bp is a necessary condition. 

We conclude that the necessity part of Theorem 1C.II is true, that is, if problem (jC.ip has 
a solution u G E^’^ with data {f,g,uo), then (IC.3P and ()C.4p are satished. We next prepare 
the proof of the existence part. 

C.2. Interior regularity and initial conditions. From [DHP031 Theorem 8.2] we deduce 
that for B G {D,N} and E G BT there exists g,B > 0 such that the realization 

fj-B — ^B with domain 'D{Ab) = {u € E) : = 0} in Lp{Q.-, E), 


has maximal regularity of type Lp(R+; E)). Thus the operator 

k-B + dt — Ab ■ okFp (R+; Lp{f^', E)) n Lp(M_|_; P(Ab)) — >■ Lp(]R+; Lp{Q‘, E)) 
is invertible for B G {D,N}. Here and in the following, the space oIFp(R+;T') is the closure 


of C'^(]R+; E) in Wp (M+; E). For A:<s — l/p<A: + l with k G No, they consist precisely of 
those functions with vanishing initial traces dfu{0) =0 for 0 < j < k, see [Ama09l Theorem 
4.7.1]. By using |Dor93[ Theorem 2.4] and a spectral theoretic argument as in |MW11] we 
may even allow for G {—Xq,oo), G (0, oo). 
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In order to obtain higher regularity results we consider the spaces 

:= - AB)-^Lp{n;E), \\u\\^, := \\{fiB - Ab)^u\\l^^^,e) for k € Nq. 

These spaces can be easily characterized by 

X^ = {u£ E) : = 0 for 0 < j < A: - 1}. 

By commuting the operator ^b + dt — Ab with {^b — AbY it follows that ^b — Ab has 

maximal regularity of type Lp(M+;Xg) for every B G {D,N}, k € Nq, that is, 

fiB + dt-A: o 1T;(M+;X|) nLp(M+;X|+^) ^ Lp(M+;X|) 

is a topological linear isomorphism. Moreover, the map 

e + df. oWl,+\R+-E)^oWl,{R+-E) 

is a topological linear isomorphism for every e > 0 and every / G No, see e. g. |MS12j . Hence, 
by commuting fiB + dt — Ab with e dt, we obtain the following result. 

Lemma C.5. Let fiD £ £ (0, oo), I G Nq, k G Nq, B G {D,N}. Then the map 

fiB + dt-AB-. oWI,+\R+-,X^e) n oWj,{R+, X|+1) ^ obL^(M+; X^) 

is a topological linear isomorphism. 

We next comment on function spaces for the initial data. From |Ama09l Section 4.9] we 
derive the characterization 

(Xix|+i)i_i/p,p = {u€ iyp2^+2-2/P(0; E) : = 0 for 0 < j < A: - jB/2 - 3/2p}. 

Then the temporal trace operator 

7t: iy;(M+;X|) nLp(M+;X^+i) ^ {X^b, X^b"\-i/p,p 
is a bounded and surjective and therefore 

+ AB,7t) : lT;(M+;X^)nLp(M+;X|+^) ^Lp(M+;X|) x (X^, X|+^)i_i/p,p 
is also a topological linear isomorphism for B G {D, X}, k G Nq. 

C.3. Boundary conditions. We will use the following result for constructing a function 
with prescribed boundary conditions (1C.51) . 

Lemma C.6. Let I G No, A: G N, p G (1, oo), let 7^ := (d^-)lr in the sense of traces and let 
X T) := Lp(r; E)) n oWj,(R+-, E)) for m G N. 

Then 7^: oE^’^(IK+ x H) —>■ x T) is a retraction and the operator 

Bi,k := ( 7 ° = 7D,7^ = lN,ll ... ,7f ■'): oE'’'=(M+ x F!) ^ TT • n x T) 

is a retraction. 

Proof. In the case H = M”, I = 0 we infer from |Aman9] Theorem 4.11.6] that 

_ 2/c_1 

Bo,k ■■= X M") ^ ll, ^ x R-i) 

is a retraction. Let Bq ^ denote a co-retraction for Ho,fc. 
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In the case Q = M” , / G No we use the fact that (e + dty : F) (I^+; F) is 

invertible for every e > 0, j G N, s G [0, oo) and every Banach space F of class T-LT. Therefore 
a co-retraction is given by = (e + dt)~^BQ ^(e -|- dtY. 

For bounded smooth domains we define such operators by a localization technique. It is 
well-known (see e. g. |GT01l Section 14.6], |PS13j l that the tubular neighborhood map 

X'. [x^t) ^ X + tvY{x), r X (—ii, R) —7> Br{T) := {x G M” : dist(3:,F) < ii} 

is a homeomorphism for some R > 0. Let {I7j : j G /} be a finite open covering of F in M” and 
let {(fj : j G /} C C'^(F) be a partition of unity subordinate to {Uj D F : j G /}. Then there 
exists r G (0, i?) such that Br{T) is covered by {Uj : j G /}. For given y G C'“((—r, r)) with 
0 < X < 1 and x(t) = 1 for |f| < r/2 we extend ipj to R"' by means of ipj{X{x,t)) := ipj{x)x{t) 
for (x, t) G F X (—r, r) so that supp ipj C Uj and d]{^ipj = 0 near F for all m > 1. 

In addition, let Uj = Bj.{x^^y with x^^'> G F for some r G (0, R) and choose rigid transfor¬ 
mations : X !->■ x^^^ + QjX with Qj orthogonal such that Qj{—en) = There exist 

iOj G C^(R"'“^) with u}j{0) = |Va;j(0)| = 0 such that for 9j{x',Xn) ■= (x',x„ + u}j{x')) we 
have Uj nQ = Uj D Sj(0j(R")) and thus Uj (IT = Uj D Hj(0j(Fo)) with Fq := R"’“^ x {0}. 
Let us construct smooth diffeomorphisms Qj of R” such that Uj H Tl = Uj PI ©^(R”) and 
UjHT = Uj n 0j(R"'“^ X {0}). Given r G (0, R/2), 'ip G Cp°{B 2 r{Q)) with ip = 1 on Bp(0), let 

^ . f y{x) [Sj(6lj(x',0)) - XnUr(Sj(6»j(x',0)))] (1 - ip{x))'Bj{x) for |x| < 2r, 

Hj(x) for |x| > 2r. 

If r G (0, R/2) is sufficiently small, then Qj is a diffeomorphism since dxQj{x) —>■ Qj as r —>■ 0, 
uniformly on R"'. Moreover, Qj has the asserted properties and satisfies —dnQj{x',(P) = 
0)) and (9™0j(x',O) = 0 for all m > 2 and x' G Bp(0). 

Choose smooth cut-off functions ipj G Cy^{Q~^{Uj)) with ipj = 1 on 0J^(supp(/9j) and 
define the multiplication operator Mj : u ipjU. With the pull-back 0* : u uo Qj and the 
push-forward 0j *: u uo 0T^ we define a co-retraction for Bi^k by 

_ 2k,—\ 

Kk9 ■■= for g G x F). 

By means of the chain rule. Holder’s inequality and the mixed derivative embeddings, it 
can be shown that the linear operators g i-a <fjg, 0^, Mj and 0j* act continuously in the 
relevant spaces and the properties of Qj and ipj with respect to the normal direction imply 
that indeed — 9- This concludes the proof of Lemma fC.61 □ 

C.4. Proof of Theorem IC.ll We have already discussed the necessity of the regularity 
conditions and the compatibility conditions on {f,g,uo). It remains to prove the uniqueness 
and existence of a solution u G for given data {f,g,uo) subject to these conditions. 

In order to prove uniqueness, it suffices to consider the most general case / = 0, A; = 1, 
where E^’^ = W/{J;Lp{TL)) n Lp{J-,Wp{Tt)) and (/, 5 ,uq) = 0. If further /tq is sufficiently 
large, then the general result of |DHP03| implies that fiQ — A has maximal regularity of type 
Lp(R+; Lp{Tl)) and this yields tt = 0 in case g,B > /^o- 

Next, we employ spectral theory to cover the case g,B G (—A^, 00), where Xq = Ao(—Ag) > 
0 denotes the smallest eigenvalue of —A^. It is well known that, since T>{Ab) is compactly 
embedded into Lp{Tl), the spectrum of A^ is discrete and consists only of eigenvalues with 
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finite multiplicity. As for Lemma lA.41 we see that its eigenfunctions belong to 
P(A^) = {ue : 7i3(A^)|r = 0 on L for j < m - l} 

for every m G N and hence belong to Wp{^) and an integration by parts implies that the 
spectrum of Ag is contained in (—oo,—A^]. A result of Dore |Dor93l Theorem 2.4] implies 
that fiB — Ab has maximal regularity of type Lp(M+; Lp(f2)) for each € (—A^,oo) and 
this ensures uniqueness. 

Existence. We construct a solution u = v} + ^ such that 

{dt + fJ-B- =■ /\ = Uh 

{dt +HB - =■ = 9j - = 0, 

{dt + fiB- A)u3 = f-f-f\ 7^A^«3 = 0, 7>3 = 0, 

for all i, j with 0<i</ + A: — 1 and 0 < j < A: — 1. This means that we first construct 
G E^’^ with prescribed initial data Ui. Then we construct with prescribed boundary 
data gj — and we finally we construct with prescribed interior data f — f'^ — f^- 

Here the functions Ui and gj are dehned according to (IC.4ah and (IC.5p by 

(C.12) Ui := f + {A - fiB)ui-i for f G N n [1, Z +/c - 1], 

(C.13) gj :=-'yBA^~^f + {dt +fJ.B)gj-i ioic j eNr\[l,k - 1], go ■■= g. 

Construction of . Let be a common co-retraction for the restriction vq: VLp(M"') —)■ 
VLp(n) for all t G [0,2A:] (cf. [AFn3t Theorem 5.22]). With the co-retraction 5""^ for the 
operator ..., 7^^^”^) from Lemma lB.31 we dehne 

:= rnS^{r^uo, r^m ,..., 0) 

+ ruS^-^iO ,..., 0, rf^ui,rf^ui+k-i)- 

Here we consider the identity operator I: 'D{I) X with X = 'D{I) = Wp^{W^) so that 
the first summand of belongs to {t i-)- e“*)H[/C'°°(M+; H4p^(f2)) E^’^. In the sec¬ 

ond summand, we consider the operator I — A: T>[A) —>■ A in A = Lp(M"') with domain 
I?(A) = VLp (M"') so that r^ui G V/^{k — l/p,p) and thus ... , 0, ..., r'^ui^k-i) £ 

VLp+^(M+; Lp(M"'))nLp(M+; Therefore belongs to E^’^, satishes '^lu^ = Ui and 

depends continuously on the data (/, uq) with respect to the norms induced by the regularity 
conditions on (/, uq) and the compatibility condition (IC.4ap . 

Construction of . From ()C.6p it follows that 'jlgj — 'jBA^Ui = 0 for all i,j G No with 
j<k — 1, i<l + k — j — 1. Thus gj — 'JbA^u^ belongs to 7solE^’^~A Near F we can split the 
Laplacian into A = Ar -|- HyBu + with the Laplace-Beltrami operator Ar = divr Vr and 
some Hr G C'°°(r). The operator Ar commutes with since it only depends on tangential 
derivatives. Therefore the normal traces hj := (j G {0,... , 2A: — 1}) of the desired 

solution G qE^’^ are uniquely determined by requiring that h 2 j+jg+i = 0 and 

7B(Ar -I- Hx'dy + = gj — for 0 < j < A: — 1. 

With the co-retraction from Lemma IC.bI we define u'^ := Bfj^{ho ,..., h 2 k-i) which depends 
linearly and continuously on {f,g,uo) with respect to the norms induced by (IC.3p and (IC.4p . 
Construction of . Finally, we shall construct a function 

G oW^+"(M+; Lp(L!)) n oIH^(M+; P(A^)), 
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which solves the equation {dt + fiB — ■= f — , where := {dt + ^B — A)u”^ 

for m G {1,2}. The compatibility conditions yield 'Jb^^ f = 0 for 0 < j < k — 2 and thus 

Hence = (9t + — A)“^/^ is well-dehned by Lemma IC.5I and the map i-a is 

continuous. The proof of Theorem 1C.II is complete. 

Proof of Corollary 1C.31 With the same arguments as above and Theorem IA.7I we see that 
fiN — Aat has maximal regularity of type Lp(M+; Lp^o(^))- Hence problem (jC.ip has at most 
one solution within the space Eq^. Analogously as for Lemma IC.51 we conclude that 

dt + m-^N-. oWl,+^{R+; Lp^om n oTH^(M+; P(A^)) 

is a topological linear isomorphism. For the proof of existence, we modify the above construc¬ 
tion of the solution n = -|- u^. For i G {1,2} we may replace u® by u® — n® G Eq^, 

since = 0 and ri®(0) = uq = 0. Then we obtain P{t) = f{t) + |F||n|“^5(t) = 0. Hence 
:= {dt + UN — A)“^/^ is well-dehned in Eg^. □ 
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